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We present an analysis of the behaviour at late-times of linear field perturbations of a 
Schwarzschild black hole space-time. In particular, we give explicit analytic expressions for the 
held perturbations (for a specihc Amultipole) of general spin up to the hrst four orders at late 
times. These expressions are valid at arbitrary radius and include, apart from the well-known 
power-law tail decay at leading order (~ a new logarithmic behaviour at third leading 

order (~ Inf). We obtain these late-time results by developing the so-called MST formalism 

and by expanding the various MST Fourier-mode quantities for small frequency. While we give 
explicit expansions up to the hrst four leading orders (for small-frequency for the Fourier modes, for 
late-time for the held perturbation), we give a prescription for obtaining expressions to arbitrary 
order within a ‘perturbative regime’. 


I. INTRODUCTION 

The study of linear field perturbations of black hole background space-times is important for several purposes. 
For example, for the investigation of the (linear) stability of black holes, the effect that black holes have on helds 
propagating in their neighbourhood or the binary inspiral of a black hole and another compact, astrophysical object. In 
1955 and subsequent years the equations describing linear field perturbations of a non-rotating (Schwarzschild) 
black hole space-time were decoupled and separated. This rendered the equations treatable semi-analytically as the 
full perturbation may be obtained as a sum of Fourier modes, with the radial part satisfying the so-called Regge- 
Wheeler (RW) ordinary differential equation. It was not until 1972 that a similar feat was achieved by Teukolsky [311] 
in the case of a rotating (Kerr) black hole space-time. In the Schwarzschild limit, the radial part of the Teukolsky 
equation reduces to the so-called Bardeen-Press-Teukolsky (BPT) equation (TU], which was also obtained in 1972. 

The RW and the BPT radial equations are satisfied by different quantities (different combinations of field com¬ 
ponents and their derivatives) and their solutions have been studied thoroughly, both numerically as well as with 
asymptotic analyses. Of particular interest for this paper, is the result by Price |3 S] (obtained by studying the field 
perturbations without Fourier-decomposing) for the behaviour at late-times of a RW field perturbation of any spin 
of a Schwarzschild black hole. Price found that its radiative Umultipole decays to leading-order in the form of a 
power-law: where t the Schwarzschild time. The analysis in [? ] of the Teukolsky equation in Kerr shows that 

the same leading-order power-law tail behaviour is satisfied by the radiative multipoles of BPT field perturbations 
in Schwarzschild. In this paper we present analytic expressions for the behaviour of general-spin field perturbations 
in Schwarzschild up to the first four orders for late-times, revealing a new logarithmic behaviour in the third leading 
order (~ t“^^“^lnt). Furthermore, although most analyses (though see, e.g., [? ] for an exception) that give the 
radius-dependent coefficient of the leading-order power-law have been constrained to large radius, our results for all 
four orders are valid for arbitrary radius. We obtain this late-time behaviour both for the field quantities satisying 
the RW as well as the BPT equation. We obtain these results by developing a method which is valid not only at 
late-times, but it may be used to obtain results valid in principle at any time regime. We now briefly introduce this 
method. 

In 1986, Leaver m derived various analytical representations for the solutions of the RW and BPT equations in 
terms of infinite series of special functions. A series of Japanese researchers later ‘revamped’ some of Leaver’s series 
representations and derived other new series representations for these radial solutions [T2HI4] . These latter series 
representations, to which we shall refer as MST expansions, are naturally adapted to carrying out small-frequency 
expansions. Small-frequency expansions yield the late-time behaviour of the full linear perturbation after integration 
over frequency. We note, however, that the MST series in principle converge for any value of the frequency, although 
the speed of their convergence decreases as the magnitude of the frequency increases. 
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The MST method is a powerful method which only relatively recently researchers have been starting to use in order 
to obtain results in black hole perturbation theory. For example, in the case of a spherically-symmetric background, 
the MST formalism has found applications in the calculation of the self-force [TS] on a point particle [TMT5] . post- 
Newtonian coefficients and gauge-invariant quantities and dynamical tidal interaction of compact objects |24) . 

The MST method has proven particularly useful for calculating the retarded Green function (GF) of the wave equation 
satisfied by the field perturbation, which is a fundamental quantity as it determines the evolution in time of any given 
initial data. The Fourier modes of the GF posess poles (the so-called quasinormal modes) and a branch cut (BG) in the 
complex-frequency plane [5^. It is known (e.g., [53120]) that the small-frequency part of the BG is the contribution to 
the GF that gives rise to the late-time behaviour; the non-small-frequency part of the BG contributes to the behaviour 
at earlier times (see [HETHao]). 

In this paper we derive in detail some of the results on the GF Fourier modes and black hole perturbations that 
we briefly presented in the Letter |3Ij . Namely, we derive a small-frequency expansion of the MST series in general, 
and of the BG in particular, which we then use to obtain the late-time behaviour of the GF and field perturbations. 
We obtain the late-time behaviour up to the hrst four leading orders at arbitrary radius in Schwarzschild space-time 
and for general-integer-spin of the field. In [T5] we used the calculation of the late-time GF derived in m in order 
to find its contribution to the self-force on a scalar charge in Schwarzschild. The MST method has also been used to 
calculate the quasinormal mode contribution to the GF in Schwarzschild space-time m and in Kerr space-time [32] : 
in the former case the GF calculation was applied to obtain the scalar self-force and, in the latter case, to obtain the 
radiation emitted given a specific perturbation source. The calculations of the quasinormal mode series in mi 131] 
involved evaluating the MST series at frequencies with ‘arbitrarily ’ large magnitude. 

Apart from an explicit small-frequency/late-time analysis, we shall also generally develop the MST formalism 
in Schwarzschild space-time. In particular, we shall present MST expressions (valid for general frequency) for the 
solutions of the RW equation for general integer spin. To the best of our knowledge, these expressions are new for 
spin-1 since the MST formalism has not yet been presented for the RW equation for spin-1 (it only has been for 
spins-0 and -2). Furthermore, we develop, also for the first time in the literature, the MST formalism for calculating 
the contribution to the GF of Fourier modes along the BG. We also give the relationships between the BPT quantities 
and the RW quantities via the so-called Ghandrasekhar transformation. We then give explicit expressions for the main 
MST quantities up to the first four leading orders for small frequency, for general spin and multipole number i. We 
note that alhough small-frequency expansions have been given for some of these quantities already, that has typically 
been done within a post-Newtonian framework, and so for the radial solutions expanded about radial infinity. In here, 
instead, we give small-frequency expansions for the radial solutions which are valid for arbitrary radius. 

For the reader who is not interested in the details and who just wants to use our results in order to obtain the 
late-time behaviour of the Green function or a field perturbation to high order, the main result is in Sec VIIA In 
particular, in Eq.(|7.4|) 


we give the late-time behaviour of the Gmodes of the GF of the RW equation; above that 
equation we indicate how to find the coefficients appearing in the equation; below that equation we indicate how to 
obtain a similar late-time expansion in the BPT case. Of course, if one wants the late-time behaviour of some given 
initial data for the field perturbation, one should convolute the obtained small-frequency expansions for the GF with 
the initial data as in Eq.(7.6). 

The layout of this paper is as follows. In SecjTTjwe present the RW and BPT equations, expressions for their GFs 
and the BG contributions, as well as the relationships between the RW and BPT quantities. In Sec |III| we develop 
the MST formalism both for the RW and BPT equations, including our new derivation in the specific case of the RW 
equation for spin-1. In Secs jl V| and [Vj we give explicit expansions of the various MST quantities, except for the radial 
functions, up to the first four leading orders for small frequency. In Sec |VI| we extract the small-frequency expansions 
for the radial functions using a novel method. In Sec VII we calculate the late-time behaviour of the GF (Sec jVII A ) 
and of a perturbation response (Sec VIIB) and compare it with highly-accurate numerical results. In AppfAj we plot 
the small-frequency expansions of MST quantities and check that they match with an independent method (Resented 
in m) which is valid in a ‘mid-frequency’ regime. In App|^we relate the radial coefficients of the solutions of the 
RW and BPT equations. 


In this paper we use geometrized units: c = G = 1. We shall use a bar over a quantity to indicate that the quantity 
has been made dimensionless via an appropriate factor of ‘2M’ (except where otherwise indicated), where M is the 
black hole mass; e.g., w = 2Mu> indicates a dimensionless frequency w and f = r/{2M) a dimensionless Schwarzschild 
radius r. 


^ We do not give the expansions explicitly for polar gravitational (Zerilli) perturbations nor for positive BPT spin but these can be derived 
directly from, respectively, axial (RW) gravitational perturbations or negative BPT spin, both of which we do derive explicitly. 










3 


II. BRANCH CUT GREEN FUNCTION 


A. Regge-Wheeler equation 


Decoupled and separated equations for linear field perturbations of a Schwarzschild black hole space-time were 
derived for axial - also called ‘odd’ - gravitational perturbations (spin s = 2) in [T], for electromagnetic perturbations 
(s = 1) in [2II], and for scalar perturbations (s = 0) in Oil]. All these integer-spin-field perturbation equations can 
be written compactly as one single partial differential equation, which in Schwarzschild coordinates reads: 




,2M 


df, = Ss 


( 2 . 1 ) 


where Xg = Xg(t, r, 9, (j)) is a scalar function that describes the field perturbation of spin s created by a source Sg and 
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( 2 . 2 ) 


is the Klein-Gordon operator in Schwarzschild space-time, where A = r(r — 2M) and M is the mass of the black hole. 
Separating variables we may obtain a complete set of solutions of the form 




(r, Cc) 


(2.3) 


where (. and m are, respectively, the multipolar and azimuthal numbers and Yi^{9,(j)) are the scalar (spin-weight 0) 
spherical harmonics. As in [33) . we treat Eq.(2.1) as a scalar wave equation, although in the electromagnetic and 


gravitational cases we have to be aware that the non-radiative {i < |s|) modes would drop out when constructing the 
electromagnetic and gravitational potentials and so these modes would have to be included separately. The radial 
functions satisfy the ordinary differential equation 


/^-d \ _2A-d uj^r 
r2 dr V dr / r^ dr A 


£{£+!) 2M(I-s2) 

n T Q 


S - g S^ryi , 


(2.4) 


where gSim = sSt.m{r,uj) are the correspon ding modes of the source Sg. We shall refer to Eq.(|2.4[ ) as the (radial) 
Regge-Wheeler (RW) equation and to Eq.(2.1) as the 4-dime nsion al RW equation (as per [33] 1. Introducing the 
standard ‘tortoise’ coordinate r* = r -|- 2M ln(r/(2M) — l), Eq.(2.4) may also be written as 


dr* 2 


- -A [<!.{£+1) Y 


2M(I-s2)' 


s^im — nsSini- 


(2.5) 


The case of polar - or ‘even’ - gravitational perturbations was derived in [3 El; the corresponding radial equation is 
the so-called Zerilli equation. Solutions of the Zerilli equation can be obtained as linear combinations of solutions and 
their radial derivatives of the RW equation for s = 2 |M] . 


We define the retarded Green function (GF) of the 4-dimensional RW Eq.(2.1) as the solution of 


V“Vc 


,2M 


sGret(a^, x') = —S(x, x') =- -5{t — t')6(r — Y)S(fl — ft'), 

r r 


( 2 . 6 ) 


that obeys appropriate causal boundary conditions, where x and x' are two space-time points. Here, H is the solid 
angle of the 2-sphere. For notational simplicity, we use time translation invariance to henceforth take t' = 0. We may 
then use the symmetries of the space-time to write 


1 


ooo-\-ic 


duj 


,Gret(x,x') = ^ / ^ Yim{0A)yL{(^'.^') e-*-‘.G,(r,/;a;), 

rr J-oo+ic 


(2.7) 


where c > 0, and the Fourier modes gGg{r,r';uj) must satisfy 


dr* 2 


2 A / 2M(1-, 

- A-f-^- 

\ r 


Gi{r,r'\uj) = -^5{r - r') = -5{r* - r*), 


( 2 . 8 ) 
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wher^ A = + 1). Correspondingly, 


sGi{r,r'-uj) = - 


fi{r<,uj)gi{ry,w) 

W{uj) 


(2.9) 


where r> = max(r, r'), = min(r, r'), and the Wronskian is given by 

W (uj) = W[fe{r,uj),ge{r,uj)] = fe^ 


( 2 . 10 ) 


Here, the ‘ingo ing’ fg = f£{r,uj) and ‘upgoing’ gi = gi(r,ui) radial functions are solutions of the homogeneous version 
of the RW Eq.(2.5) which, for general spin, behave asymptotically as 


h 




( 2 . 11 ) 


where r* = r*/(2M), and are complex-valued constant coefficients (we give higher-order terms of Eq.(2.11) 
in Appj^, and 


It is then easy to see that 


r*-)■- 1 - 00 . 


W = 2iujAi". 


( 2 . 12 ) 


(2.13) 


It will be convenient for the next section to define similar ingoing and upgoing solutions of Eq.(2.5) without choosing 
a specific overall normalization: 










, -OO, 

-I- f* -)► -1-00, 


(2.14) 


and 


X7 






Y-up,tra 

e-ZV 0 O , 




CX), 

+ 00 , 


(2.15) 


where ^j^up,mc/ref/tra ingoing/reflection/transmission coefficients. It is clear that 

\^in _ ^in,tra n v'up _ -^up.tra 

s^£ — s^£ Jli s-^£ — 9^' 


Using the standard spherical harmonic addition theorem, we may now rewrite Eq. (2.7) as 

OO 

sGret{x,x') = -^^^^'^{2£+l)Pe{cos'y)sGY’\r,r';t), 

where 


e=o 


-1 poo-\-ic 

,G^°*(r,r';f) = 7 ^ dw sG£{r,r ;oj)e~ 

J — oo+ic 


(2.16) 


(2.17) 


(2.18) 


and 7 is the angle between the spacetime points x and x'. 

The upgoing radial solution g£, unlike fi, possesses a branch cut (BC) in the complex-frequency plane starting at 
the origin a; = 0 and extending down the negative imaginary axis [TTJ [531 inn . This BC is inherited by A™ and by the 
Fourier modes sG^(r, r';w) of the GF. This BC in gg, however, only occurs as a change of sign in its imaginary part 
as the frequency crosses the negative imaginary axis, therefore \gg\ and |y4“| do not possess a BC. We will use a new 


^ We note that our definition of A is slightly different from that in El 

^ Note that due to a typographical error the right hand side of Eq.(2.3) of Ref. m should read — f'^gi instead of gif'^ — fig'^ so that 
it is indeed W = as claimed below Eq.2.7 |27| . 
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‘frequency variable’ cr > 0, so that when oj = —ia then it lies on the negative imaginary axis and when oj = +ia then 
it lies on the positive imaginary axi^ We define 6A{a) = Aj^{—ia) — A-^—iu) for any function A = A(uj) possessing 
a BC along the NIA, where A±[—i^ = limg_>o+ A(±e — ia), with ct > 0, where a = 2Ma and Co = 2Moj. 

The contribution from the BC to the f-mode is given by 

OO 

sGf^{r, r'-1) = ^ J da SsGi{r, r'; a)e~'^*. (2.19) 

0 

Using the obvious symmetry gi(r,—uj*) = g|(r,w) together with the boundary condition gi{r,-\-ia) ^ as 

f* —>■ OO and the fact that all three functions gi+{r, —ia), gi-{r, —ia) and ge{r, +ia) satisfy the same homogenous linear 
second-order differential equation (namely, the RW equation) with real-valued coefficients along both the negative 
and the positive imaginary axes (since lj = ±ia with cr > 0 there), it follows that 

Sgi = iq{a)gi{r, +ia), (2.20) 

for some real-valued function q{a), and that W[gi-{r,—ia),gij^{r,—ia)] = —2iaq{a). We refer to q{a) as the ‘BC 
strength’. The symmetry fg{r, —oj*) = f^{r,uj) together with the fact that has no BC means that fg{r, —ia), with 
cr > 0, is a real-valued function. Putting all these results together, we find that the discontinuity along the BC of the 
Umodes of the GF is given by 


6sGe{r,r'-,a) = -2ia^^fi{r,-ia)fe{r',-ia), 
\W\ 

which is a purely-imaginary quantity. The reader may refer to [281135] for details. 


( 2 . 21 ) 


B. Bardeen-Press-Teukolsky equation 


The Newman-Penrose formalism offers an alternative way of describing spin-field perturbations of a Schwarzschild 
black hole background space-time to that provided by the RW formalism (one key advantage of the Newman-Penrose 
formalism, however, is that it generalizes to Kerr space-time). The scaled Newman-Penrose scalars = 4';s(t, r, 0, (p) 
obey the following equation [HHin]: 


V“V„ 


2s 

„2 


-r + 


Mr'^\ d 


d 


COS0 d 

A J dt ' '' ' dr ' “ sin^ 0 dip 


+ (r - M) — +i 


= T„ 


( 2 . 22 ) 


1 — s cot^ 0 
2 

where Tg is the matter source term. The scalings of the Newman-Penrose scalars here are given by 

{ ipo or s = -1-2 or s = —2, 

(po or r'^(p 2 , s = -1-1 or s = -1, 
ip, s = 0, 

where '0o/4 are the radiative Weyl scalar s, (po / 2 are the radiative Maxwell scalars and (/j is a massless scalar field. We 
note that in the scalar case (s = 0), Eq.(2.22) is the same as the 4-dimensional RW Eq.(2.1). 

As for the RW equation Eq.( |2.3[ ), we may obtain a complete set of solutions of Eq.(2.22) in the form 

'^sit,r,0,(p) = e““*sW„(0,(^)si?£m(?’,a;), 

where aY(.m{0,4') are the spin-weighted spherical harmonics 

Yi^{ 0 ,(P), 


(2.23) 

[37]. These are given explicitly by oh£m(0,())) = 


-_iYi^{0,cP) = 


±2Wm(^;^) — 


1 


vW+l) ■ 


Tde - -^d^ 
smO 


Yim{0,(p), 


^{i-me+ !){£ + 2 ) 


di — cot t 
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2i 

sind 


T 


2i cos 0 
sin^ 0 


> 1 , 


W^(6»,(;i), e> 2 . 

Sin 0 


(2.24) 


The spin-weighted spherical harmonics satisfy the following relations: 


^ We note that in |27l 1281 i31| (and in the BC literature references therein) we used a different symbol for the frequency a = iu. The 
symbol used there coincided with the symbol for the ‘renormalized angular momentum’ parameter introduced later on that is used 
throughout the MST literature and which we also use in this paper; hence the reason for the change of symbol to a. 
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• the conjugation relation sYl^{9,(j)) = (-l)'"+'*_sy'z(_r„)(6'»; 

• the orthonormality relation 




i' ^mm' 1 


• the completeness relation 


OO ^ 


Y. Y sYl^{e,ct^)sYi^{9',cl>') = 5{cose-cose')5{ct^-cl>')-, 

^^|s| m——l 


• the parity relation sYim{T^ - 9 ,tt + (j>) = (-l)^_sYf(_m) (0, 4>)\ 

• and the generalised addition relation 


YsM0,^)s.YU9' 


(2.25) 


where 


a = tan ^ (cos 6* sin 0' — sin d cos 9' cos[(j) — (j)'), — sindsin(^ — (/)')), 
(3 = tan“^ (sin0 cos 0' — cos 0 sin 0' cos{(j) — (j)'), — sin0' — cj)')). 


In particular, when si = S 2 = s, 


YsYU0,^)sYL{9',^') = ^(-1)^ (^^|^)''pi!f^(cos7)e-(“+«, (2.26) 

■m ' 


where Pj,°'’^\x) are the Jacobi polynomials. 

The radial part of the functions in Eq.(|2.23|) satisfies the ordinary differential equation 


A—+2(r-M)(s + l)—+ 
dr^ dr 


d f {r'^UJ — 2is{r — M)) 


+ Aisujr + s(s + 1) — A 


sRim = r^sTem, (2.27) 


where sTim = sTimir^w) are the corresponding mode s of t he source Tg. We shall refer to Eq.(2.27) as the (radial) 
Bardeen-Press-Teukolsky (BPT) equation and to Eq.(2.22) as the 4-dimensional BPT equation. We may write the 
former in self-adjoint form as 


(ir'^UJ-2is{r - M))r'^uj ' 

A®+^ + A® R --h 4zswr -b s(s -f 1) - A 

dr dr \ A 


sRim = r^A® sTtm, 


or, writing sRem(r,uj) = (A /r)gXemir,uj), 


'' d^ 

dr 2 


n 2is(r - 3M)uj , A 2MA 
■OJ H-^s- - -A-r - 


/\s/2+l 

s\im — sRdm- 


(2.28) 


(2.29) 


Le t us denote by sXi ^ general homogeneous solution of Eq. (2.291. Asymptotically, the homogeneous version of 
Eq. (2.291 takes the form 


r j 2 

I / ■ 

—^ -b (W — ISK) 

dry 


''d^ 
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-bw^ -b 
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sXi ^ 0, 

sXi ~ o> 


r* —>■ —OO, 
f* —)■ -boo. 
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where k = 1/(4M) and correspondingly the solutions behave (omitting dimensionful constant factors) as linear 
combinations of 




/ A 


\{2My 


±s/2 


±2u;r* 


— 00 , 


+ 00 . 


We may define ‘ingoing’ solutions of the homogeneous version of Eq. (2.29) by their asymptotic behaviour as 


in 

sXi 


-s/2 

V2M 7 

j^tn,inc^s ^—iujr^ _j_ Sg+2a;r* 


(+-!)■ e 


sJ-i-i 


r» —>■ —oo, 

f* -)► +00, 


(2.30) 


and the ‘upgoing’ solutions as 


sX7 


j?" 


'■{2M) 


l+s ("77 _ (2M)^”® ("77 _ ^ 


j^up,tra^-s^+iuir. 


\2M 


\2M 


(2.31) 




r* —)■ + 00 , 


for general BPT spin. The corresponding ‘ingoing’ and ‘upgoing’ solutions of the homogeneous version of the BPT 
Eq.(2.27) respectively behave asymptotically as 


Til 

glXp 


TjV, 

Tjin,inc _i _ 

s^e 


A ®e r* —)■ — 00 , 

^-l^-iuir, ^J^n,ref ^-i- 2 s^+iuir, ^ + 00 , 


(2.32) 


and 


j^up 

s^i 


sRg 


up,inc 4-ia;r, 


p' 

sR-e 


up,tra i- 2 s„+i 


sRe 


•np,ref ^-s^-iuir. 


— OO, 

+ 00 , 


(2.33) 


where g^pg ^^jjg incidence/reflection/transmission coefficients of the ingoing radial BPT solution; sim¬ 
ilarly for fgp ^]pg upgoing solution. It is convenient to define the following ‘ingoing’ and ‘upgoing’ 

solutions and coefficients with a hat on, which are normalized with respect to the corresponding transmission coeffi¬ 
cients: sRT — sRT/ sRT’*^‘^’ sRT’ — sR^^/ bRi^’*^°', Then it is easy to see that 


W^{oj) = [sW,sR7\ = A^+i [sR^tsRT' - sRTsRf) = 


(2.34) 


where primes denote differentiation wrt r. The quantity W'^ is a ‘generalized Wronskian’ in the sense that it is 
independent of r. Similarly to the RW case, the solution sR^^ has a BC along the negative imaginary axis of the 
complex frequency plane, which is inherited by IT^, whereas si?7 BCs m- 

The GF of the 4-dimensional BPT equation is the solution of Eq.(2.22| with the source Tg replaced by the distri¬ 
bution 6{r — r')6{cos6 — cos 9')S{(j) — 7)/-\77(- H may be expressed as 


gGret,T{x, x) = ^ A7r')sG™‘’^(r, r; t)sYtm{0, ip)sYl^{9', tp'), 

i,m 


(2.35) 


where 


GT^^ir,r';t) = 


f'OC-\-ic 


du} sGJ (r, r'; uj)e 


sGj{r, r';uj) = - 


,+(r<,a;),i?7(r>,a;) 


' —oo+ic 


w^ico) 


In order to obtain an expression for the BC contribution gGf to the £-mode of the GF of the 4-diniensional 

BPT equation we proceed similarly to the previous subsection for the RW equation. This contribution can be expressed 


^ We keep a spin s subindex in the solutions and of the BPT equation while we did not for the solutions and of the RW 

equation because of the explicit spin-dependence in the asymptotic Eqs.( |2.3^ and | |2.33| l for the former set of solutions (as opposed to 
Eqs.l |2.lT| l and ( |2.12[ | for the latter set); this leads to the explicit spin-dependence in Eq.l |2.3^ . 



















as in Eq.(2.21) for the RW case, but with sGf^ and 5sGi replaced by sGf"’^ and SgGj, respectively. We then note 
that = sRg^ (^)^) and that A“®_siZ“^(r,+ i(t) goes like jr as f —>• oo (neglecting a constant 

factor). It can be shown that all three functions sRi+{r,—icr), sRt-{r,—ia) and +iCT) satisfy the same 

homogenous linear second-order differential equation (namely, the BPT equation) with real-valued coefficients along 
the imaginary-w axis (since w = ±ia with cr > 0 there). The two former solutions have the same large-f behaviour, 
which is linearly independent from that of the latter. Therefore, it follows that 


SsRg^ = iq'^{a)A 


"-sR'^^{r,+ia), 


(2.36) 


for some real-valued function q^{(j), and that sRi^, 




= —2iaq^(a). We emphasize that the ‘BC strength’ 


function is calculated using the ‘upgoing’ radial function along the positive i magin ary axis with the opposite spin- 
sign, as opposed to the same spin-sign for calculating q in the RW case (see Eq.(2.20)). The property sRi^{r, —oj*) = 
sR™ (f, together with the fact that has no BC means that sRp{r, —iu), with cr > 0, is a real-valued function. 
Putting all these results and properties together, we find that the discontinuity along the BC of the f-modes of the 
4-dimensional BPT GF is given by 


5sGj (r, r'; cr) = -2ia'^-^sRTir, -ia)sRT{r', -ia). (2.37) 

\WT\ 

In the next subsection we relate the BPT quantities to the RW quantities and, in particular, we present an alternative 
way (namely via the RW ‘BC strength’ q{a)) of calculating the BPT discontinuity 6sGj. 


C. Relationship between RW and BPT quantities 


The so-called Chandrasekhar transformation relates, in the homogeneous case, solutions of the RW equation to 
solutions of the BPT equation. We note that while the RW Eqs.(2.1) and (2.41 are symmetric under s O —s, the 


BPT Eqs.( 2.22) and (2.271 are not. We here write the Chandrasekhar transformation compactly for spin s = 0,-1, —2 
(see [3H] for spin-2 and, e.g. [33] for spin-1; there are similar transformations in the case of positive spin - which of 
course only changes the BPT equation, not the RW equation - but we do not deal with these in this paper). 

Let us generically denote by sXe{t,r) a homogeneous solution of the 4-D RW Eq.(2.1) after factorizing out the 
angle-dependence via scalar spherical harmonics; similarly, we generically denote by s^e(tj f) a homogeneous solution 
of the 4-D BPT Eq.(2.22) after factorizing out the angle-dependence via spin-weighted spherical harmonics. The 


corresponding Chandrasekhar transformation is then 




dr - 

r ^ t 


Id 


dd-i 


(2.38) 


up to a normalization constant. We already gave in Eq.(2.24) the angular counterpart of the above transformation. 


i.e., the transformation from the angular factor in the 4-D RW Gmodes (namely, the scalar spherical harmonics) to the 
angular factor in the 4-D BPT Gmodes (namely, the spin-weighted spherical harmonics). We note that if Eq.(2.24) 
were naively applied to the modes £ < |s| it would yield the zero function. 

It is useful to write the Chandrasekhar transformation in the frequency domain explicitly for each spin s = 0, — 1 
and —2 separately. We now generically denote by and gRe homogeneous solutions to the (radial) RW Eq.(2.4) 
and (radial) BPT Eq.(2.27), respectively. Introducing the operators 


"4 


f{r)’ 


= ^ _ 

dr /(r)’ 


where /(r) = Ajr'^, we may express the BPT solutions in terms of the RW solutions as 


oRe = -oXi, 

r 

-iRe = ci(r - 2 M)V_iXe = dr [/(r)_iX; + tw.iXi ], 

_2Rt = C2{r-2MfV^{r_2Xi) 

= C2 [ 2 (zr^a; -I- r - 3M)r(/(r)_2A( -|- zw_2A^) -|- {£{£ + l)r - 6M)/(r)_2W] , 


(2.39) 
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where cq, ci and C2 are constants of proportionality and primes denote differentiation with respect to r. Conversely, 
we have 


o^£ — — 
Co 


(2.40) 


_,.Y, = ^ 

Aci V r 

-2^e = 


1 

Aci 




fir) 

2 f —iRl 


((£ - l)e(e + 1 )(£ + 2 ) - 12 iMw)c 2 
1 

((£ - !)£(£ +!)(£ +2)- 12iMuj)c2 




2(3M — r — iujr^) 
r-2M 


-2Re+ 


(£(£ + 1) - 2r^uj^ + Sicor + 4) - 2Mr (£{£ + 1) + 9ia;r + 10) + 24M2 


-2^1 


r{r — 2 M )2 

We will now use the Chandrasekhar transformation in order to relate the BPT and RW Wronskians as well as 


the BPT and RW ‘BC strengths’. The specific normalizations (2.11) and (2.12) of the RW radial functions fi and 


ge yield specific normalizations for the corresponding BPT radial functions via the transformation Eq(2.38). We 


will denote the BPT radial functions and coefficients with these specific normalizations with a tilde superscript, i.e., 
si?™ = Oi^(r)fi, = Oai{r)ge and 5 ??™’**^“, si?“^’*’'“, .. .are, respectively, the coefficients si?™’*’'“, ... of 

si?™ and ■ In Appll^we use the Chandrasekhar transformation in order to find these BPT radial coefficients and 
Wronskian in terms of the RW radial coefficients and Wronskian. We are here using the obvious notation of O^j (r) 
for 0{t,r) after replacing ‘dt’ by ‘—iuj’ in it. 


Let us now try to find an alternative e xpress ion to Eq.(2.37) for the BPT 6 sGj in terms of the RW q{cr). We 
first re-express the BPT modes gGj of Eq.(2.35) in terms of 
transformation Eq.(2.38): 


sGj{r,r';a) = -Oujir<)Ouj{r>) 


the RW solutions fg and gg by using the Chandrasekhar 
feir<,oj)ge{r>,u}) 


j^in.tra f^up.traTjrn-' ' 
s -tt£ s ^ 


Using this expression we can find the discontinuity of the BPT Umodes as 


dsGj{r,r';a) = 


W 




0^(r)Ouj(r')SsG£(r,r';a) = - 


2iaq{a)sR, 


in^tra 




(2.41) 


«i?r(uw)si?r(r',a;). (2.42) 


In deriving Eq.(2.42) we have used the fact that W/{sR!' iP'' ■ sRf^’ does not have a BC (as can be seen from 

the expressions in App|B| ) and we have made use of Eq.(2.21 1 . 

Comparing Eqs.(2.37) and (2.42) we immediately obtain a relationship between the RW and BPT ‘BC strengths’: 


q^{a) = 




in,tra 


TD‘ 

fyup^tra 

S-tt£ 


q(a). 


(2.43) 


Using the results in App|^ we have that 


Q^i^) = <li^) ‘ 


2 cra; ’ 

/3^(a++ 2M2/3+) 


{2a^a* 


■2M2/3*)’ 


s = 0 , 
s = - 1 , 

s = -2. 


(2.44) 


and 


^in.tra 


w 


1 

4M2 


|^T|2 |W|2 |W|2 


2fJOioo 

(a+ + 2M2/3+)2 

2MV2/3oo 


s = 0 , 
s = - 1 , 

s = - 2 . 


(2.45) 
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Inverting the relationship and using the results in App|^ we explicitly obtain: 


g 

|iy|2 


= c 


9 


T 




with 


(2.46) 


C = 4M^ 


C = -- 


IQ\M^ 


(l - 4cr + 4(7^) , 


C = 


(A - s2 + 1) 

2bQM^(\-2aY{l-af 

A^ — 2A — 6a 


s = 0, 
s = -1, 

s = -2. 


(2.47) 


The small-frequency behaviour of C for s = 0 and —1 is already manifest in the exact result above; for s = — 2 we 
expand it as 


C = 


256M6 

(A-2)A 


6(A2-2A-l)d (l3A4-52A3-hl6A2-h72A-h36)dM ^ 

1- + - - - — -fO(d3), s = -2. 


(A-2)A 


(A - 2)2A2 


(2.48) 


From Eqs.(2.37), (2.46) and the fa ct that the radial functions and are generically of the same leading order 
(order zero) in d as d —)■ 0 (see Secs IF t^I), it follows that the RW and the BPT BC modes SsGe and dsGj are of 
the same leading order in d as d —>■ 0. As a consequence, the RW and BPT GFs are of the same leading order in t as 
t —>■ 00 , as we explicitly see in Sec VII That is, the RW and the BPT quantities describing black hole perturbations 
decay at the same rate at late times. 


III. MST FORMALISM FOR THE RW AND BPT EQUATIONS FOR GENERAL SPIN 


The MST method for the Teukolsky equation in Kerr (and, therefore, for the BPT equation in Schwarzschild) was 
given in ST for general spin, and earlier in [T^] (henceforth MSTa) just for spin-2. To the best of our knowledge, 
however, the MST method for the RW equation has only been given for spin-2, which was done in |13j (henceforth 
MSTb). Therefore, the MST method for the RW equation still has not been developed for spin-1 (obviously, the RW 
spin-0 case is essentially just the same as in ST with s = 0). In this section, we develop the MST method for the RW 
equation for general spin: for spin-2 we recover MSTb, for spin-0 we essentially recover ST with s = 0 and, for spin-1, 
to the best of our knowledge, the results are new. We also develop the MST method for the BPT equation which, 
although already existing in the literature, will allow us to emphasise the connections between the MST formalism for 
the RW and BPT equations for general spin s = 0, —1 and —2. In particular, we write the expansions for the radial 
solutions of both the RW and BPT equations in terms of just one set of ‘universal’ coefficients (namely, a„). When 
referring to the literature, we shall use the generic term MST to refer to all MSTa, MSTb and ST. 

We shall we use the notation of X, 


homogeneous RW solutions ^'^d the BPT solutions gR 

of the RW and BPT equations when using the specific normalization as in MST (i.e., the one in Eqs.(3.12) and (3.20) 


and (i.e., with a slight change in the subindices with respect to the 


respectively) for the ingoing/upgoing solutions 


below). A similar change in the subindices notation applies to their incidence, reflection and transmission coefficients. 


A. Series of hypergeometric functions 

In this section we shall assume that s < 0. For the RW functions we write 

X'Y = e“f"+ie-*‘‘^’'*p“(a;) = (f- l)-“#"(a;) = (1 - (-a;)-“ p“(a:), (3.1) 

where a; = 1 — f. This leads to the ordinary differential equation 

x{l — x)py'{x) + (l — 2iw(l — xY — 2(1 -I- s)x)pY'{x) (3.2) 

-I- ((£ - s)(£ -I- s -I- 1) -I- 2id;(l -|- s)(l - x))pY(x) = 0. 

Correspondingly, for the BPT functions we write (as in ST) 

rY = (f- l)-"e-“"*p“(a;) = (f- l)-"-“p“(x) = (-cc)”"”"*p“(a;). 
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which leads to the equation 

a:(l — x)p™''{x) + {AiMujx{l — a;) + 2{2iMuj — l)x + (1 — s — AiMu)))p™' {x) (3.3) 

+ {(.{i + 1) + 4iMw(l — s)(l — x)')p™{x) = 0. 

In terms of p™{x) and p™{x) the Chandrasekhar transformations take the form 

pII'(x) = CoPo^(x), 

= cip“/(x), 

P-2(^) = C2P-2"(^), 


and conversely 

pil'(x) = -pIP(x), 

Co 

= £(£ + l)ci + (1 - 2*w(l - x)^) p“i(x)] , 

p‘^ 2 (^) = 77 - 77 - 77 - 7 - 7 — [ 2:(1 — x) ( 2 x + 1 — 2 ia;(l — x)^} p%'(^) 

— {£{£ + l)a;(l — x) + 4w^(l — x^ + 6iu>(x + 1)(1 — x)'^ — 2{2x + 1 )) p™ 2 (*)] ■ 
We now follow MST and introduce the expansion in terms of hypergeometric functions, 


S"W = ivr E a /*-" - - + * ty +»+1 - X 

n——(x> 


p“(x) = ivr £ 


r(i - 2 iuj) 

2^1 (—n — ly + s — iuj, n + v + s + l — iui; 1 — 2iuj; x), 
r(— n — V — iui)T(n + + 1 — iui) 

r(i - s - 2 itd) 

2 Fi{—n — V — iui, n + v + 1 — iui] 1 — s — 2iui] x), 


(3.4a) 


(3.4b) 


where is a normalization constant which we specify later on. Here, the parameter v is referred to as the renor¬ 
malized angular momentum and is determined by the requirement that these series converge both as n —>■ oo and 
n —>■ —oo, it has the property that either lu = £ + 0{ui) or i/ = —£ — 1 + 0(ui); see MSTa or ST for a full discussion. 
These series representations in terms of hypergeometric functions converge Vr € [2M, oo). 

In these terms, the Chandrasekhar transformations follow from the standard identity (Eq.I5.5.1 [ID]) 

2.^1 (a, b] c]x) = — 2 Fi{a -I- 1 , 6 -|- 1 ; c-|- l]x), 
da; c 

for some parameters a, b and c. Inserting into their respective differential equations and using the standard hyperge¬ 
ometric function identities 


X 2Fi{a, b] c; a:) = 


a{b — c) 


{a — b){a — b + 1) 
c{a -I- 6 — I) — 2ab 
(a — b — l)(a — & -I- I) 


2 ^ 1 ( 0 -I- l,b- l;c;a;)-|- 

b{a — c) 


2 Fi(a,b;c;x) + 


(a — b — l)(a — b) 


2 Fi(a -l,b+ l]c;x), 


x(l-x)^2Fi(a,b;c;x) = '2Fi{a + l,b- l]C]x)+ 

da; (o — o)(a — 0 -I-1) 


a 6 ( 2 c —a — 6 — 1 ) , ab(a — c) , 

2Fi{a,b]C]x) + 7-TTT- 1 --y 2^1(0 -l,b+ l]C]x), 


(a — 5-I-l)(a — & — I) ’ ’ ’ (a —5)(a —5—1) 
we find that d„ and a„ must satisfy the same three-term recurrence relation: 


Q^n^^n+l “t“ 'yn^n—1 — 0; 


(3.5) 
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where 


iuj{n + V + 1 — iu}){n + + 1 + s — iuj){n + + 1 + s + iui){n + v) 

{^Ti -\- 2z/ -t- 3) 

= -A(n + v){n + t/ + 1) + ((n + v){n + + 1) + 

iijj(n + V + ibj){n + v — s + iuj)(n + v — s — iuj)(n -\- v + 1) 


In = 


(2n + 2 j^ — 1) 


The overall normalisation of the coefficients in the homogeneous Eq.(3.5) is, of course, irrelevant for the value of 
the a„ but the above form has the advantage that all denominators are bounded away from 0 in the perturbative 
(small-frequency) regime. By ‘perturbative regime’ we essentially mean the frequency regime where v is real - see the 
end of Sec |IV| for further details. We choose the normalization do = oq = 1, so then we have d„ = a„, Vn S Z, and 
from now on we will write down all series using the ‘universal’ set of coefficients a„. These coefficients a„ are equal, 
for s = —2, to the aj) in MSTb as long as the same normalization is chosen for the two sets. 

We could alternatively choose to include the T-functions in the coefficients, that is, write 


Rw T(—n — ly + s — iu})T(n v + s — iio + 1) 
a„ oc Qn - 


T 

Ciri 


r(i - 2 iQ) 

T(—n — V — iw)r(n + v — iuj + 1) 


(3.6) 

(3.7) 


r(i - s - 2 iuj) 

where the constants of proportionality are independent of n and so just reflect the normalisation of the series. A 
particularly convenient choice is 


,RW _ 


= a„(-jz -I- s - iuj)-niv + s - iu} + 1)„, 


a„ = ani-iy - iuj)-n{v - iu + 1)„ = a, 


(v — iu + 1), 
’ (ly + iu + 1), 


(- 1 )". 


(3.8) 

(3.9) 


As we choose the normalization gq = 1, we also have aj = = 1. Using this convention the corresponding 

three-term recurrence relations have coefficients 




T 

= - 


n + iy + 1 — s + iu 

i Td I ~o;n, 

n + iy + l + s — lu 

n + ly + 1 + iu 

i Td ~o;n, 

n + ly + 1 — lu 


= /3n, 

Pu = Pn, 


In 


7n = 


n + ly + s — lu 

I i ~"inj 

n + ly — s + lu 
n + v — iu 


n + ly + lu 


—'In- 


(3.10) 

(3.11) 


Note that, in the perturbative regime, where iz is real, Eq. (3.111 may be reexpressed as 

B'^ = 0* 

h'n t-’ni 


qI = at 


T * 
'In = 'In- 


Up to irrelevant overall normalisation, the coefficients a'^, 0'^ and yj and corresponding are the same as the 
corresponding quantities in Eq.l23 ST. 

As the corresponding coefficients differ by a scaling that tends to 1 for large |n|, ly is the same for RW and for BPT. In 
particular, v depends only on |s| which can be seen directly since under the transformation s —>■ —s, n+v —>■ —n—iy — l, 
0n is invariant while a„ and 7 „ are simply interchanged. 

As the event horizon of the Schwarzschild black hole is approached, we have 


Xl’ 


-p“(0) 


and 




(r- 1) 


—icur* in 


P“(0), 


2M. 


That is, our solutions (3.4al and ( 3.4b[ ) are normalised according to (see Eqs.(2.14) and (2.32|) 

r(— n — V + s — iu)T{n -I- ^ -I- s -I- 1 — iu) 


dr) 


Xl^’tra ^ e*‘^p”(0) = ^ 

n— — oo 

oo 

j^n,tra ^ (2M)2«p“ (0) = (2M)2« ^ 


r(i - 2iu) 

r(— u — V — iu)T(n V 1 — iu) 


r(l — s — 2iu) 


(3.12a) 

(3.12b) 


We note that the particular normalization choice, 

N'U = 


r(l — s — 2iu) 


r(—iy — ia;)r(l -i- ly — iu) ’ 

yields the specific normalization used in ST for the ingoing BPT solutions, and so that is our choice henceforth. 


(3.13) 
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B. Series of Coulomb wave functions 


An alternative expansion is useful for the construction of the ‘up’ solutions. In terms of the variable z = uir = ujr, 
the RW equation may be written as 

f ^ Mz) = 0 . 


z — UJ 

— id) 


Writing sXt{z) = (l — 7 ) hs{z) this becomes 


z^h" + \z^ + 2ijjz — l{l + 1 )] hg = oj 


z{h'g + hs) — (1 — 2ioj)h!g — (s^ — (1 — -hg — ujhg 


(3.14) 


The left hand side is the operator defining the Coulomb wave equation (Eq.33.14.1 [ID]) with solution satisfying ‘up’ 
boundary conditions at infinity given by 

H+i-Ld,z) = W,^^),+ i{-2iz) = e^^{-2izY+^U{£+l-iCd,2i + 2,-2iz), 

where z) denotes the (unnormalised) irregular Coulomb function (Sec.33.2(iii) [JD]) and W denotes the Whit¬ 

taker function (Sec. 13.14 [ffl), We use a ‘hat’ to denote that in writing the above we have dropped the conventional 
normalisation prefactor where ai{—u}) is the Coulomb phase shift, which is irrelevant to our cur¬ 

rent discussion. U denotes the irregular confluent hypergeometric function (Sec.13.2 [ID], 'k in the notation of [H]). 
Again, following Leaver El and MST, this suggests that we introduce the expansion for the upgoing solution 


h^p(z) = iv;p £ 


r(n -I- 1 / -I- 1 -I- s — iuj)T{n + v + \ — ioj) - ^ 


r(n -I- 1 / -I- 1 — s -I- iuj)r{n + v + 1 + ioj) 


yTa„iL++^(-w,z) 


= N^Pe^-{-2iz)’'+^ 


CLn, X 


n— — oo 


r(?T, V -\-1 + s — iuj)T{n v -\-l — ioj) 
r(u-|-:/-|-l — s-f zw)f (n -I- ^ -I- 1 -I- iu)) 

(—2iz)”[/ (n -I- -I- 1 — iw, 2n + 2v + 2, —2iz). 


(3.15) 

where N'^p is a normalisation constant that we will specify later so that our normalisation agrees with ST. Inserting 
into Eq.(3.14|, eliminating second derivatives using the differential equation satisfied by IL^(— 77 , 2 :), and noting the 
following identities which follow from standard properties of the functions U : 


(—V, z) H-- zHt(—ri, z) + i—.—iHf , (— 77 , z), 

{L+1){2L + 1) ^ L{L + 1) ' L{2L + 1) 


1 fr+( \ ■ (L + l-iv) 

-HZ{-ri,z) = -i 

— H^i-n z] = i + ^- ’iy) ■^+ /_ N 

dz (L + 1 )( 2 L + 1 ) ^+ 1 ^ 


T+^ ^ I + ’^V) ■£r + 


L{L + 1) 




L{2L + 1) 




we find that a„ must satisfy the same three-term recurrence relation Eq.(3.5). 
Similarly, the BPT equation becomes 


sR'i{z) -I- (1 -I- s) 


’ +1) .g,(») + (1 + _ '’(< + 1) -»(»+ P \ = 0 . 


Z — UJ z 

— S — lUJ 


{z-^f 


r(z — Lo) 


Writing gRi{z) = z ^ * (l ~ 7) hg{z) this becomes 


z^h'g [z^ -I- 2 ( 0 ; -I- is)z — i{l -I- 1)] hg = oj 


z{h'g + hg) — (1 — s — 2iu})h'g — (1 — iuj){s - 1 - 1 - iuj)-hg + i{s + iui)hg 


(3.16) 


The appropriate Coulomb function is now (—w — is, z) and the corresponding expansion for the upgoing solution 
is 


= (-i)*iv-p £ 


r(n -I- 7^ -I- 1 -I- s — iuj)T{n + v + 1 — iuj) - 


r(n -I- 7^ -I- 1 — s -I- 7a;)r(77 + v + \ -\- iui) 


- 'is,z) 


= (-l)*lV;Pe-(-2zz)'^+i ^ ^ 

n——oo 


r(7T, -I- 7 ^ -I- 1 -|- s — iuj)T{n -I- 7 ^ -I- 1 — iuj) 


dn, X 


77-1-77-1-1 — S-|- 7w)r(77 -|- 77 -|- 1 -|- ILo) 

(—27z)"17 (77 -I- 77 -I- 1 -I- s — tw, 277 -I- 277 -I- 2, —2iz). 
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The Chandrasekhar transformations in this case follow term by term from the Whittaker function identities: 

^l+ii5,L+i 


('-!) 


(‘to) 

—iuj \ 

1 W^TO.L+i(-2*^)j 

i^TTo 

1 

1 —I 

I 

II 

(‘to) 

—iGj \ 

1 W^TO.L+i(-2*^)] 



2—iuj 


W2+i<5,i+l(-2iz), 


or equivalently 


^ ~ x) ~ 2) = - (1 - ^) H+^A-^ + hz) 


(1 - ^) (z (1 - ^) iJ++^(-w,z) j = z (1 - H++^{-uj + 2i,z), 




where 


2 ^0—w( ^+*1^- 




-1 


Using the relationship between a„ and aj we may write our solution in the alternate form 

h"P = (-l)"ie-""e*"(-2z)"+i y (^ + 1 + ^ - al{2iz)^U (n + ^ + 1 + s - 2n + 2t^ + 2, -2iz) , (3.17) 

_ + 1 - S + iuj)n 


where we have made the choice 


iV^P = U 


_ 1 r(t^ + 1 - s + iu})T{v + 1 + w) 


r(i/ + 1 + s — iuj)T(iy + 1 — ioj) 


so as to agree with the normalisation of ST. However our original form serves to highlight the boundary conditions 
and the link to the RW solution. 

The corresponding upgoing BPT solution is given by 


Ry = 


(/ % \ S Z(jJ 

1 -^) .TO. 


(3.18) 


This series representation in terms of irregular confluent hypergeometric functions converges Vr > 2M. Since 
U(a,c,x) ~ x~°‘ as \x\ to 00 , it is straightforward to write down the asymptotic forms 

^up ^ ^—l—2s^i{z-\-oj\n.z) 

where 


HC = (2*)-^e-5^e-*5('^+i)2-i+“ V (^ + 1 + ^ - + 1 - 

^ + 1 — s + + 1 + 


(3.19) 


Clearly, then, from Eqs.(2.33) 


s Av A(jj In (jj 


^up,tra ^ (^2iyA'ie 
j^up^tra ^—l — 2s^Qhiuj 


(3.20) 


To conclude this subsection, we note that we could have used the ansatz sXt{z) = (l — ks{z) for the RW 
equation, giving 


z^k" + [z^ + 2 a;z — + 1 )] kg = to 


z{k” + kg) — (1 + 2iuj)k'g — (s^ — (1 + iw)^) —kg — uik. 


(3.21) 
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and the upgoing solution 


k7{z)= Y. 


n— — oo 


r(n + — s + 1 — iuj)r{n v + 1 + iuj) 
r(n + + s + 1 + iu))T{n + v + 1 — iCj) 




= e*"{-2iz)"+i Y 


r(n + — s + 1 — iuj)T{n + v + \ iuj) 


r(n + + s + 1 + i(jj)T(n + z/ + 1 — iuj) 


■:r 7 an{—‘ 2 iz)'^U (n + z/ + 1 — zw, 2n + 2z/ + 2, —2iz). 


While this expansion seems to more naturally capture the boundary conditions of the RW equation, the Chandrasekhar 
transformation is less natural in terms of it, so we will discuss it no further. 


C. Relation between the two solutions 


Using the standard relation Eq.15.8.3 [ID], we may reexpress in terms that are better suited for discussing its 
behaviour at radial infinity: 


= iVr(l - x) 


iujx 


Y 


r(—n — 1 / + s — iuj)r{2n + 2z^ + 1) 
r(n + z^ + 1 — s — z(D) 


a„.x 


(1 — x)’^ 2 ^i ( —n — zz + s — iuj, —n — v — s — iuj] —2n — 2u] 


I — X 


+ Nril-x) 


— v+iuj ^iujx /jj.'j— 


r(n + z^ + 1 + s — ztLi)r(— 2n — 2z^ — 1) 

/ . pTT X 

^ 1 (—n — V — s — IUJ) 


(1 — x) ^ 2 Fi (n + z^ + l + s — iuj, n + z/+l — s — iuj; 2n + 2v + 2; 


1 — X 


The second term can be obtained from the first by the substitution n —> —n, v —> —v — 1 and correspondingly the 
terms are denoted by Xq and respectively, so 


X in I V —^1 

s ~ ^0 -^0 


(3.22) 


(Note that is invariant under this transformation.) It is easily checked that each of these terms is independently 
a solution of the RW equation and moreover are linearly independent (MSTb). In terms of z = ujf = a)(l — x), 


, OO 

( 2 \ —lU) / 2. \ + _ 

5-D (s) H 


n——oo 


In an identical fashion we can write 


with 


r(—n — z/ + s — iuj)T{2n + 2i> + 1) 
r(n -\- V + 1 — s — iuj) 

2 F 1 n — V — iuj, —n — v — s — iux, —2n — 2v, —^ . 


r(—n — V — za;)r(2n + 2z^ + 1 ) 

_Li -^ 

l(n + z^ + l — s — lUj) 


—n — V — IUJ, —n — V — s — luo] —2n — 2z/; — 

z 


a„x 



Din 

iLg 

= Ro 

("--1 

\ / 2 ; ' 

. u+iu) 

Vo; 

/ Vw. 

) 


(3.23) 

(3.24) 


ilTM 


(3.25) 


This series representation for i?™ converges Vr > 2M. 

To obtain solutions in terms of confluent hypergeometric functions suitable for discussing the behaviour at the 
horizon, we introduce the auxiliary solutions involving the regular Coulomb wave function 

Ti(-w, z) = 1 (-2iz) = +l-ii^,2£ + 2, -2iz) 

= + 1 + z^, + 2, 2zz), 
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where we include the prefactor for reasons that will become clear when we consider the Chandrasekhar transformation 
below. The function Fg{—i],z) then satisfies the following identities 

1 A / N (L + 1 — iri) A . ^ r? A . ^ {L -\- irt) ~ 

-Fr,{-n.z) = ~ h(2t + 


d ^ ^ ■ L{L+1-iri) r] 

S = -•(L + i)(2L + i)'^‘+‘<-’''"> + WTY) 


c 2 \ •(d' +1)(^ + 'i'n) ^ ^ 

Fl{-V,z) - , .V — FL-i{-'n,z). 


L{2L + 1) 

Proceeding as before we construct the corresponding solutions: 

r(n + + 1 + s — iui)T{n + ly + 1 — iut) 


X‘, 


( ! ) \ — 

1--) (-2*zr+i^ 


zJ ' ' —' r(n + + 1 — s + ia;)r( 2 n + + 2 ) 

n— — oo ^ ^ ' 

( 2 iz)”M (n + ly + 1 + iui, 2n + 2iy + 2, 2iz). 


(3.26) 


and 


RUz) = (l-^) 


x —s—iu) 




r(n + i^ + l + s — iLu)T(ji + ly + 1 — iui) 
r(7z T “t“ 1 “t“ ?ca)r(2/z T 2iy -t- 2) 


where 


{2iz)^M {n + V + 1 — s + iuj,2n + 2iy + 2, 2iz ), 

^ r(t/ + 1 - s + iw)T{iy + 1 + w) 

2r(z/+ 1 + s — ia;)r(j^ + 1 — iw) ’ 


(3.27) 


(3.28) 


for agreement with Eq.l39 ST when z > 0 . Our normalization for Xq follows from that of Rq via the Chandrasekhar 
transformation, and so for s = —2 it does not coincide with the normalization choice in MSTb. Specifically, for 
s = —2, our Xq is equal to that in MSTb times 

T{iy + 1 — s + iuj)T{iy + 1 + iQ) 

T{iy — 1 — ioj)T{v + 1 — iui) 

The Chandrasekhar transformation when using these series representations follows term by term from the Whittaker 
function identities: 

1 ~ “J ^0 n 1 — — j ^iuj,L+l (—2iz) j = (L + 1 + zw) — —J 1 (—2iz), 


2 




(l-^) M^^j^^i{-2iz)\ = {L + l + iuy){L + 2 + iuj)z{l-^'^ ^ 2 +io,L+i(-2*z), 


or equivalently 


((^~ z) ^n+A-^,z)j=(l-‘^'^ F+^^{-u; + i,z), 


1 — iuj 




( / • —tuj \ / / J \ ^ 

1 --) = z[i--) F+^^{-Q + 2i,z), 


where 


Fo — uj \ — + j ^1 — —^ 


^dz V z- 

These solutions can be related to X^p and R^^ using the identity in Eq.6.7(7) Vol.l [12] (valid for b ^ Z): 


e“"C(o, b, 2iz) + - a, b, -2iz), Re(z) > 0, 


M(a, 6 , 2 zz) = 


m 

r(6 -a)^ 


t/(a,5,2iz) + ^e-(“-^)"e2*"C(&-a,&,-2iz), Re(z) < 0, 
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the first and second terms on each line yielding the incoming and outgoing wave solutions at infinity, respectively. 
Thus, X^{z) = X!l{z) + X'iiz), where 


/ lD\ 

Xl{z) = (l - - j {2izY+^ ^ 


r(n + ^ + 1 + s — luj) 

;-rrX 

r(n + ^ + 1 — s + lUj) 


(—2iz)"17 (n + z/ + 1 + iuj, 2n + 21^ + 2, 2iz) 


(3.29) 


— _ 

X'^iz) = (l - {2izY+^ ^ 


n— — oo 


r(n + 1 / + 1 + s — ia;)r(n + z/ + 1 — iui) 
‘r(7T, + zy + 1 — s + *w)r(n + z/ + 1 + i(D) 


(—2zz)"17 (n + z/ + 1 — zw, 2rz + 2zz + 2, —2iz ), 
where the signs correspond to Similarly, Rq{z) = R\{z) + R''_{z), where 


EYY) = (l - 


— S — IU) 


-^-Y2izY+^ f2 


{—2izYU {n + V + 1 — s + iuj^2n + 2v + 2, 2iz ), 

V —s—iQ ^ 

1--J z-^-Y2izY~^^ J2 


r(n + + 1 + iuj) 

(3.30) 

r(n + z/ + 1 + s — iuj)r{n + z^ + 1 — zw) 


n= —oo 


‘r(7z + zz + l — s + za;)r(n + zz + 1 + zo;) 


{—2izYU (n + zz + l + s — zw, 2n + 2z/ + 2, —2iz). 


Note that our naming convention here follows ST and is opposite in sign to MSTb. In particular, the minus solutions 
are just multiples of the corresponding ‘up’ solutions. 

Critically as noted by MST, The functions X^ and Xq solve the same differential equation and have the same 
analytical behaviour as functions of Z] similarly for i?Q and Rq. Therefore, they must be proportional: 


X^ = K^XY and R’'^ = K^R^ 


(3.31) 


where K^, is the constant of proportionality. Equating the corresponding Laurent series we can obtain explicit 
expressions for in terms of v and our a„ coefficients. The results are given for general spin by ST (based on their 
Teukolsky equation analysis), we repeat them here for completeness specialised to Schwarzschild space-time: 


K, = 


e“(2d;)'’-'^-’'"’- 


2-Sjrr(i g 2zd;)r(r + 2i^ + 2) 


r(r -I- zz -I- 1 — s -I- zw)r(r -|- zz -|- 1 -|- za;)r(r -|- zz -|- 1 -|- s -|- iuj) 


E 


E 


\n— — oo 


r(7z -I- r -I- 2zz -I- 1) r(rz -|- zz -|- 1 -|- s -|- iuj) r(zz -|- 1 -|- iw) 

(zz — r)! r(rz -\- v 1 — s — iCj) r(zz -|- 1 — zw) 

1 (zz -I- 1 -I- s — ioj)n (y + \ — iLo)r, 


(r — rz)!(r + 2^ + 2)„ (zz -|- 1 — s -|- zu;)„ (zz -|- 1 -|- iuj)r 


(3.32) 


The parameter r here is an arbitrary integer number. 

The above expressions lay out the foundations for taking the f —)■ oo limit of the ingoing solutions and thus finding 
their incidence and reflection coefficients. First we relate the quantities at —zz — 1 to those at zz. From Eq.(3.28) we 
have 




sin ( 7 r(zz -h zw)) sin ( 7 r(zz + iuj - s)) 

^ 7 *”" Y" •> ?" C C 1 \ ^. 

sm ( 7 r(zz — tu})) sm ( 7 r(zz — lui + s)) 


Then, from Eqs.(3.29) it follows that 


XT 


sin (7r(zz -I-zw)) _j 
sin (7r(zz — iw)) 


(3.33) 


(3.34) 


and 




(3.35) 
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where the up per/l ower s ign c orresponds to Re( 2 :) positive/negative. 

From Eqs.(3.22) and (3.29) and 13.7.3 [40] we can take the z —>■ oo limit and obtain the incidence and reflection 


coeflicients of the ingoing RW solution: 




and 


where 


A'; = 


sin (7r(iy — zw)) 

/2 i^(iy+i) o-i-iiz r(iz + 1 — s + iiD)r{iy + w + 1) r(n + 1/— zw + 1 + s) 


n— — oo 


r(rz + v + iuj + 1 — s) 


(3.36) 


(3.37) 


(3.38) 


We can proceed similarly for the BPT solution. By taking z —)■ c» in Eq.(3.24) we obtain the incidence and reflection 
coefficients of the ingoing BPT solution: 




and 


where 


j^in^ref ^ — l — 2s 


sin (7r(zz + s — iu})) 

( jy' I • \ Ai^ Ind) 

(Ay + ze A_y_ijA_e , 




^(^+i-s) 2 -i+s-io ^(‘^ + 1 - s + za;)r(z^ + zo; + 1) 
r(iy + 1 + s — za;)r(z^ — zw + 1) 




(3.39) 


(3.40) 


(3.41) 


n— — oc 


r(n + 1 / + ioj + 1) 

We note that one can obtain the radial incidence, reflection and transmission coeflicients of the RW solution from 


those of the BPT solution (or viceversa) via the Chandrasekhar transformation Eq.(2.40) (or Eq.(2.39)). In doing so, 
the leading order for large-r would be annihilated and one would require a higher order term. 


IV. LOW-FREQUENCY EXPANSION OF THE COEFFICIENTS AND OF ly 


In this Section our goal is to provide the low-frequency behaviour of the MST series renormalised angular momentum 
ly and series coefficients an- We will provide the expansions explicitly up to t he firs t fiv e leading orders. The behaviour 
of the coefficients and may b e ded uced immediately from Eqs.(3.8) and (3.9). 

We start by noting that Eqs.(3.2) and (|3.3|) reduce to the hypergeometric equation when a) = 0, indeed it was 


precisely for this reason that Leaver El and MST wrote them in this way. For the ‘in’ RW and BPT solntions we 
want the regular solutions corresponding to, respectively, 

27^1 (£ -|- 1 -|- s, —£ -|- s, 1; x) = 2 -f^i (—£ -l-s,7-|-I-|-s,l;x) 


and 


2^1 (7 -b I, -£, l-s;x)= 2Fi(-£, £+ 1 , 1 -s; x). 


The left hand sides of these expressions correspond to ly = £ and the right hand sides to = —7 — 1 when w = 0. In 
fact, under ly —>■ —v — 1, q;„ equals 7 _„ and /3n equals /3-n and therefore, a„ satisfies the same recurrence relation 
as a-n under ly —v — I; the equivalent symmetries hold for the RW counterparts (a]/'^, £3^ and y^'^) and 

for the BPT counterparts (a)^, Pn y)])). This symmetry stems from the fact that the renormalized angular 
momentum v was introduced into the Ordinary Differential Equation Eq.(I19) ST in the form iy{iy -|- 1), which is 
invariant under z/ —> — z/ — 1. In addition, we may determine the expansion of zz abou t £ from that about £ — 1’. 

With the natural ansatz that = 0(u)l”l), the 3-term recurrence relation Eq.(|3.5|) can be solved directly yielding: 


v+\- {£+\) 


A(15A - 11)-b 35-*-b 6(A - l)s2 _2 
A(4A-3)(4A+1) 




4(A - 2)A3(4A - 15)(4A - 3)3(4A -b 1) 


jLo^ + Oiui^) (4.1) 
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for the renormalized angular momentum. Note that the expansion of v must be even in s in Schwarzschild space-time 
since it also arises through the expansion of the RW equation which manifestly has this property. As it will be needed 
later on, let us define as minus the coefficient of uP' in the above expansion for ie, 


V = i — V2Up + 0{(jP'). 

For the series coefficients themselves, the above procedure yields 
(£ - s + 1)2(£ - S + 2)2(£ - S + 3)2(£ -5 + 4)2 


a4 = 


• — 3 V 

“ 24(^ + 1)(£ + 2)(2£ + 1)(2£ + 3P(2£ + 5)^^ ~ 24(£ + ip(£ + 2^(1 + 3)(2£ + 1)(2£ + 3)2(2£ + 5) 


(^-s + 1)2(£-s + 2)2 2 

02 = - 


(£-s + 1)2(^-s + 2)2 


4{£ + l){2£ + l){2£ + 3)2 4{£ + ip{£ + 2){2£ + 1){2£ + 3) 


ioj^- 




mi£ + ip{£ + 2){£ + 3)(2^ - 1)(2£ + ip{2£ + 3)4(2£ + 7) 


(£-5 + 1)2 _ 

2(£+l)(2£ + l)*‘^' 

+ 


(£ - 5 + 1)2 
2 (£+ 1 ) 2 ( 2 £+ 1 ) 


uj + O(w^), 

3 ( 12 . 6 ) 


OJ 


P'ir\£,s) 


P[f^\£,s) 


0-1 = - 


8£(£ + 1)4(£ + 2)(2£ - 1)(2£ + 1)3(2£ + 3)2(2£ + 5) 

3(12.6) 
.3 


U{£ + 1)3(£ + 2)(2£ - 1)(2£ + 1)3(2£ + 3)2(2£ + 5) 
u}'^ + O(w^), 


lUJ 


+ sf (£+ 5)2 _2 

—--W 




2 £( 2 £ + 1 ) 
+ 


2£2(2£ + l) 8(£-l)£3(£+l)(2£-3)(2£-1)2(2£+1)3(2£ + 3) 

,_,4 

8(£ - 1)£4(£ + 1)(2£ - 3)(2£ - 1)2(2£ + 1)3(2£ + 3)"’ 

2 ^^ ' "12 + 1 ) 2 (^+ S )2 _ 


lUJ 


-co^ + O(w^), 


(£ + 5-1)2(£ + 5)2 2 , _ 

4£(2£-1)2(2£+1) 4(£ - 1)£2(2£ - 1)(2£ + 1)' 


P[f!\£,s) 


+ O(w^), 


48(£ - 2)(£ - 1)£3(2£ - 5)(2£ - 1)4(2£ + 1)3(2£ + 3) 

21^-L 112^^_l „12 P_ 3 , 4 ( 15 , 8 )(£, 5 ) 


(£ + 5-2)2(£ + 5-1)2(£+s)2 ^._3 

“ 24(£ - 1)£(2£ - 3)(2£ - 1)2(2£ + 1)*^^ 

^ (£+5-3)2(£+5-2)2(£ + 5-1)2(£ + 5)2 
96(£-1)£(2£-5)(2£-3)2(2£-1)2(2£+1) 


24(£ - 2)(£ - 1)2£2(2£ - 3)(2£ - 1)2(2£ + 1) 


LO^ + O(w^), 


(4.2) 


96(£+1)(£ + 2)(2£+1)(2£ + 3)2(2£ + 5)2(2£ + 7)‘^ ^ 

- 5 + 1)2(£ - 5 + 2)2(£ -5 + 3)2 ^ ^2£ + ll) (£ - 5 + 1)2(£ - 5 + 2)2(£ -5 + 3)2 


(4.3) 


0(u 


where p(’">") (x, y) denotes a real polynomial with integer coefficients of degree to in a; and n iny and recall A = 1(1 + 1). 
The precise form of the polynomials is easily determined to high order but is too long to be useful in printed form for 
general s. For completeness, we give, as the most important example, the term in the renormalised angular momentum 
for 5 = 0, —1, —2: 

(P 0) = - 2) (3240 + 8733A - 82625A2 + 155295A3 - lOSOOOA^ + 18480A®), 

PqT^ (P 1 ) = -(^ - 2)(405 + 1917A - 3240A2 - lOdOSA^ + 19435A4 + 54495A^ - 84840A® + 18480A^), 

PqT^ (a, 4 ) = 51840 + 102816A - 953424A^ + 51522A^ + 123233A'‘ + 85775A® - 2415A® + 61320A^ - 18480A®. 

The naive pattern evident in the leading behaviour under the ansatz a„ = 0(wl"l) is 


f ( (£+1-5)„)"(£+1)„2- 
( 2 £ + 2 )„( 2 £ + l)2„n! 


naive , 
On = < 


(-iw)l”l+0(w"+4) n>0. 


((-£-s)H)'(-£)|n|2l"l 

U-2£)|„|(-2£-l)2|„||nr > ^ > 


(4.4) 
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where {x)n denotes the Pochhammer symbol (a;)„ = r(a; + n)/r(a:) = x{x + 1)... (x + n — 1). Note that reflectin g the 
symmetry noted at the beginning of the section, we ha ve o„ —>■ a-n under £—)•—£— 1. We have also derived Eq.(4.4) 
via an alternative met hod: by imposing that Eq.( |3.32' I for K,^ is independent of the parameter r. 

Critically, Eq. (4.4) reveals that the ansatz a„ = 0(a;l"l) is flawed since the denominator vanishes whenever 
\n\ >2^+1 while, in addition, the numerator vanishes when |n| > ^ + s + 1. A detailed analysis reveals that the 
correct ansatz ioi v = t + 0{uP') and s = 0 is: 


rO(a;l’^l) n>-^, 

an= -{2e+l) <n<-e, (4.5) 

n<-{2e + l). 

The same rules apply to s 7 ^ 0 with the overrides: 

s = —2 : ttn = n = —I + 1 and n = —£, (4-6) 

s = —1 : a„ = n = (4.7) 

s =+1 : a„ = 0(wl"l“^), n =—£ — 1, (4.8) 

s =+2 : a„ = 0(wl"l“^), n =—£ — 1 and n =—£ — 2. (4.9) 


Inserting this revised ansatz into Eq.( |3.5[ ), together with an expansion for u, yields equations that can be solved 
recursively along rising and falling diagonals (treating v as if it were ag) to very high order. Through the diagonal 
nature of this procedure, it becomes clear in all cases that the general terms given above work for through the 
following orders while beyond this they must be supplemented by expansions based on the corrected ansatz: 


General expressions of Eq. (4.3) valid for s < 0 to ord 


r o(w”+2^+i), 

[o(cC-"- 2 ), 


n > —£, 

-(2£+l) <n<-£, 
n < -(2£ + l). 


We may usefully turn this around, if we wish to work uniformly to order uj^, the most anomalous behaviour occurs at 
either n = —£/ — £— I or at n = —(2£ + 1) and w e m ay use the general expansions except when min(t' +1,2£—1) < N. 
For example, the expansions to order uj'^ of Eq. ( |4.3[ ) are valid for any £> 2 while we need to calculate the expansions 
for ^ < 2 using the revised, correct ansatz to order w"*: 


s = 0 : £ = 0 


£=1 


£ = 2 


ai 


ozo 

1 


1 _ 1 _ 
= -2^^+2^ 
2 7 


In 331 n 

r - 


- LUJ -r - 

360 360 


■UJ 


6 


1591 _n 31723 _o 73956143 _4 1 _ 1 

- oj H- VjJ - UJ . a _9 = —luj H- 

34020 50009400 ’ ^ 9 54 


a_i = 


—I- iuo - 

9 27 2835 


2 -2 

a-3 = ^ 0 . 


2 -3 
243 


■UJ 


1243^4 


34020 


1-4 '1 --3 

u; , as = —* 0 . 


1323 

1. 1 2 

Oi = --IUJ + -UJ - 

3 6 

1 - 1-2 
— tCJ 


13_4 3_2 1_9 r(^r 

1575 ’ ^ 50 20 67500 

19 _2 1325203 _4 

- ^ , 


a_"s =- 

® 4725 

1741 _4 
--w , 


a_i = --iuJ 

6 

25 _2 

a_-? =- UJ 

^ 722 


6 


2447 _n 2221 

- iuj H- UJ , u 

18900 18900 ’ 

169 _3 467 _4 

20520 4104 ’ 


= 1 - uj‘‘ - 1 

30 3591000 


5 .-3 

- lUJ^ 

228 


5070221 _4 

ooo/tnono^ ’ fl—4 

0.2 


“ ^ 4332 

25 .,-,3 


545 _3 5341 _4 


25992 


-ito 


3249 


-UJ 


a_ 


32840892^^ ’ “ 2166 

5 _4 5 _3 47 _4 12 _2 1 _3 

9702 ’ ^ 882 10584 ’ 245 35 

3 1 _2 21793 _3 38383 .4 „ 79 7082 

- IUJ -\ - UJ - lUJ -I- UJ , V = 2 - UJ - 

10 10 441000 1323000 210 9261( 

l._ 1 _9 1613 ._3 1399 _4 l_2^1.-3 

=-90^" + 60*" - 


Oi 


i _ i _2 IDid _ 

0_1 =- lUJ - UJ - lUJ - u 

5 10 31500 31500 

7-4 n 49 _4 

-UJ , 0-4 = 0 , 0-5 = — _ UJ . 


0-3 


"1422 


112338 


(4.10) 

9679 _4 

- UJ 

34020 
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ff ^ oo _4 iU ,_3 61 _4 Zb _2 I 

84942 ’ ^ 2079 12474 ’ 567 2f 

2 _ 1 2 351521 _3 117151 _a , 169 

ai = —lOJ H- oj - loj H- oj , V = 6 - 

7 14 12224520 9779616 630 

3 _ 1 _2 6019 _o 29887 

a_i =- luj -w- lUJ - Lo , a_2 =- oj' 

14 205800 1852200 175 

= 0 . 


'l4 

1 


a-a =- ioj^ H- 

2100 12600 


11 -4 

- -OJ , a_4 : 


2750517000 
3,13 15979 4 

'- lOj'^ - OJ^, 

70 5145000 ’ 


: 0, a_5 


s = -1 : £=l 


= 2 


'25 


1 -4 1 -3 13 -4 

a4 = , a3 = , a, 

3 3 2 14977 _3 24847 _ 4 ^ 47 .3 43908007 .4 

ai = —loj H —OJ - XU! H- OJ , 10 = ! - oj - oj , 

4 8 33600 67200 ’ 60 71064000 ’ 

17 ._3 17 _4 5 ._3 107_4 

a_i =- XU - u , a _2 =- loj - oj , 

^ 282 282 ’ ^ 47 564 ’ 

200 _2 10 _3 152936921 _4 150 .3 85 _a 48 _a 

- OJ H- lOJ - u! , a _4 = - xoj H- oj , a _5 = - oj . 

2209 141 307419903 2209 4418 2209 

47 _4 __ _ 20 _ 2 5 _. _3 440941 .4 


a-3 


ai 


5 _4 10 _3 47 _4 20 _2 5 _3 440941 

1782 ’ 411 2646 147 63 14261940 

8 8 _2 32987 ..3 56647 _a „ 169 _2 6832249 _a 

—XOJ “t“ —OJ —- XOJ “t“- OJ , 10 — 2 —- OJ —- OJ , 

15 45 330750 992250 ’ 420 74088000 ’ 

1 _ 1 _2 5707 _3 7397 _a 11 _a 

a_i =- XOJ - OJ - XOJ - OJ , a_2 

20 40 168000 336000 ’ 

7 -4 „ ■l!l -4 

“-= = - 6083 “’ “-= = "’ “-= = - 53132 “ ■ 


'6760 


OJ , 


=-2 : e = 2 


10 _4 
“^=891^"’ 


5 ,_3 47 _4 

03 = XOJ - OJ , 

^ 72 864 ’ 


15 _2 5 _3 730781 _4 

02 =- OJ - XOJ - OJ , 

49 28 6338640 


5 5 _2 12029 _3 19519 .4 

oi = —loj H- OJ - XOJ H- OJ , 

6 18 52920 158760 ’ 


„ 107_2 1695233_4 

10 = 2 - OJ - OJ , 

210 9261000 


1 --3 7 -4 

0-1 =- XOJ - OJ , 

20 40 


11 


a-2 = 


12840 


7 


a-3 = 


1926 




0-4 = 0, O-H = —- 


98 

11449' 




In order to evaluate q^{cr) via Eqs.(5.10) and (5.11) we also need the sum of the a„. Using Eq.(4.3l we obtain, to 
order oj'^, 


“ , (A + «=) 


lUJ 


(4A2 - 7A + 4) + A(4A + 9 ) 5 ^ + 6A(4A - 3)s^ + 3(2A - 1) + A - 6 ) + 2A(4A - 3)(7A - 6 )s 


P 3 


rXOJ^ + 


4A2(4A-3)2 

Pi 


OJ 


24(3 - 4A)2(A - 2)A3 48(A - 2)2A4(4A - 15)2(4A - 3)^ 


-oj''' + O(oj^), 


(4.11) 


where 

P 3 = (A - 2)A^ (4A^ + 213A2 - 309A + 108) + i + OOA^ - 81A + 54) s® 

- 18z(A - 2)A(4A - 3)5^^ - zA (36A3 - 175A2 + 432A - 252) - 12z(A - 2)A(4A - 3)(5A - 3)3^ 

- a (36A^ - 279A3 + 686A2 - 738A + 288) - 6 i(A - 2)A2(4A - 3)(7A - 6 )s, 


(4.12) 
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P 4 = -(A - 2)2a2(448A® + 13584A^ - 68220A® - 372693A^ + 2636097A‘^ - 5229081A3 + 4967622A2 - 2347380A + 437400) 

- (A - 2)A(448A® + 6288A® + 19620A^ - 334989A3 + 958311A2 - 841995A + 109350)s® 

+ 18(A - 2)2A(4A - 15)(4A - 3)^(4A2 - 102A + 45)/ 

+ (-256A^ - 24896A® + 45680A^ + 1948740A® - 12101346A’^ + 28846314A^ 

- 34305849A3 + 21625542A2 - 7333740A + 1312200)s® 

+ 6 (A - 2)A(4A - 15)(4A - 3)2(100A^ - 1194A3 + 4425A2 - 3744A + 540)s® 

+ (9600A^° - 142432A® + 532648A® + 2117966A^ - 23736786A® + 7967033^^^ 

- 130327029A^ + 113471766A^ - 54566136A2 + 14346720A - 1749600)/ 

+ 2(A - 2)A/4A - 15)(4A - 3)2(540A^ - 6334A3 + 25641A2 - 40896A + 15660)5^ 

- A(256Ai° + 57152A® - 505072A® + 7324A^ + 14370666A® - 68251182A® 

+ 151062669A'‘ - 183278700A3 + 123738192A2 - 43759440A + 6415200)5^ 

+ 2(A - 2)A2(4A - 15)(4A - 3)2(84A® - 3342A^ + 21295A3 - 47598A2 + 42048A - 11880)s, (4.13) 

with special cases 


s = 0 : £ = 0 


e = i 


e = 2 

s = -1 : e=l 

1 = 2 

s = -2 : £ = 2 


E 

n— — oo 
oo 

E 

1 = —cx 
oo 

E 

1 = —oc 
oo 

E 

1 = —oc 
oo 

E 

1 = —oc 
oo 

E 


7 7 


122 


- ILO - UJ 

9 54 945 


6343 _3 276257 


n— — oo 


n— — oo 


68040 


= 1 - 


= 1 - 


= 1 - 


1 _ 854 _2 19669 _o 

— _ UJ — _ lOJ — 


2500470 
3079726441 


9025 72200 


7982161250 


+ OiCo^), 
uj'^ + O(w^), 


1 _ 53 _2 377 

—icu - U) - 

2 882 3528 


_3 167323193 .4 , ^,. 5 , 

iuj — .^ 4 “ CJiuj ), 


4855023684 


3 _ 19643 

-VjJ H- 

4 441800 


^ 2 - 


2212413 
3534400‘ 


— 




12^^^ 3528 


59 _2 329 _3 3421065415 

UJ — _ „ IUJ — 


= 1 - 


1728 
2045 


10841072759029 
19128349520000 

uj^ + 0{uj^), 


Cd^ + Oiuj^) 


5 _ 25 _2 

-lU! - UJ - 

6 882 5292 


iuj^ — 


236995065216 
1706608417 


23750538336 


Cd^ + Oiuj^). 


Finally we should note that the power series for u cannot be valid Vo; S C. The reason is that, for example for 
w € M, the exact value of ly is known to be real for small values of w but, as w increases, u reaches some half-integer 
value and then it suddenly picks up an imaginary part (see, e.g.. Table 1 in ST; we have found a similar behaviour 
for (D on the negative imaginary axis). The power series for however, is purely real for w real and so it cannot 
reproduce this behaviour. In this paper, however, we are only interested in the small-d) behaviour, where the power 
series does converge. We have been referring to such regime as the ‘perturbative (small-frequency) regime’. 


V. LOW-FREQUENCY EXPANSION OF THE BRANCH CUT INTEGRAND 


It is clear from Eqs.(2.17l and (2.191 that the late-time asymptotics of the BC contribution to the GF is provided 
by the small-d expansion of the modes 6sGi. The small-d expansion of 6sGi, in its turn, is given by the small-d 
asymptotics of the radius-independent quantity g/|VFp and the radial solution f^. In this section we will derive the 
small-d asymptotics of g/|VFp. We will use these asymptotics later in Eq.(??) in order to prove that sGf^ ~ 
to leading order for large-t, for general integer spin. 

Specifically, in the following subsections we provide explicit expansions for small-a) up to the first three leading 
powers (which actually correspond to the first four leading orders, since the third order has the same power of w as 
the fourth order but, as we shall show, it contains a logarithm in uj) for the various perturbation quantities (except 
for the a„ and v, which we gave in the previous section, and for the radial functions, which we give in the next 
section). We provide these expansions for general integer-spin s and multipole number £. These expansions will yield 
the first three leading powers (so four leading orders) for late-times of the multipole-1! GF and field perturbations. 
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The above comments apply to the RW GF and field perturbations but a similar argument applies to the BPT ones. In 
fact, we shall give the small-(t expressions explicitly for BPT quantities; one can then readily find the corresponding 
expansions for the RW quantities via the transformations given in Sec jll C| 

From now on and for the rest of the paper we shall restrict ourselves to the case that the BPT spin is a negative- 
integer, s = 0,-1,—2 (RW spin, of course, is indistinctively positive or negative). BPT quantities for positive spin 
can be obtained from those for negative spin via the Teukolsky-Starobinskii identities [34l 1^ . 

Also, as mentioned, from now on we will focus on small-frequeny expansions up to the first three leading powers of 
the frequency. It is easy to see from the results in the previous section (e.g., compare the general-^ and -s expressions 
in Eq.(4.3) with the specific-mode expressions in Eq.(4.10)) that the expressions that we shall obtain for general 
multipole-£ and spin-s in principle are not necessarily valid, up to the first three leading powers of the frequency, for 
the three specific modes £ = 0 (and s = 0) and £ = 1 (and s = 0,-1). Indeed, some general-s and -£ expressions 
that we shall give appear to have singularities at the s and £ values for these modes. We have dealt with these three 
modes separately by carrying out small-frequency expansions after setting the corresponding values of s and £ right 
from the start. Remarkably, we have found that our general-s and -£ expansions up to the first three leading powers 
of the frequency, for and all subsequent quantities derived from these two quantities, actually give the 

correct result for two of these anomalous modes, namely for .^ = 1, with s = 0 and — 1For the mode s = £ = 0, 
the general-s and -£ expressions do not give the correct result and we present the results for this mode (as well as for 
£ = —s = 1 and 2, for completeness) in Sec VE 


A. A, 


The quantity introduced in Eq.(3.311 is ne eded in order to obtain the Wronskian below (in Eq.(5.4)). We obtain 
the following small-w expansion of from Eq.(3.32): 


Tf 2^r (£ + i)r(2£ + 2)r(i-s) __,^, 

" V^r2(£-s + i) "" 

i2^+^T{2£)T {£ + I) F(1 - s) {2{£ + l)£ {He - 2H_s +1 e)+£'^+£- s^) __ 
V^(^ + l)r2(£-s + l) 


(jj 




T{£ + l)r{£ - s + 1)"^ 


(5.1) 


where 


= 


_ 4 (6 (£2 -H f - 1) s2 -h £{£ + 1) (15£2 -h 15^ - ll) -k Ss'^) {2He_s +He- 2 H 21 - 2H2e+i + ln(2w) -k je) 


£{£ + l){2£-l){2£+l){2£ + 3) 


4 ' "^(fTl) ^ ^ ^ j ~ - 4 (Hi - 2FI-.) " - 4 j - 

2 {{£ + 1) {8£^ - 16£2 - 24^ -k 9) s"^ -f £‘^{£ -k 1)^ (8£^ -k 16^^ - - 13£ -k 4) - 


r(2) 


^2(^+1)2(2£- 1)2(2^+ 3)2 
(16€® -k 104£® -k 164£'‘ - - I3l£’^ -3£+ 18) s^) \. 


(5.2) 


Doing similarly for we obtain 

(-+2^-ir(l-s)r(£+ l)T\£ + s + l)__g+, 


A. 


-1/-1 — 


r(2£-k l)r(2^-k 2):y2 


-U) 


■0(w^+"-i) 
(- 1 )" 


(5.3) 


In obtaining Eq.(5.3), we have used Eq.(4.4) together with the fact that r(— n -k w) ~ r(ra++ as w —>■ 0 when n is a 
nonnegative integer. 


® Specifically, for, e.g., s = 0 and k = 1, the general-s and -I expression for Ku gives the wrong coefficient in the third leading power 
of u) and for it gives the wrong leading order. However, for this mode, the two incorrections in Ku and in somehow 

miraculously cancel each other out to give the right result for up to the first three leading powers of u). 
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We note that the asymptotics of in Eq.(5.1) and those of in Eq.(5.3l imply, via Eq.l 68 ST (and since 

I > |s|), that it'-i,-! is not necessary for obtaining to the first four leading orders starts playing a 

part only in t he ne xt order), except in the cases 1 = 0 and £ = 1. As mentioned at the start of this section, though, 
Eqs.(5.1) and (5.31 are in principle not valid for ^ = 0 and 1, although the general-s and -£ expressions that we give 
for quantities from now on are, somewhat surprisingly, also valid for £ = 1 (with s = 0,1). The case 5 = ^ = 0 we 


treat separately in Sec VE 1 


B. Radial coefficients and Wronskian 


We now turn to the coefficients in Eq.(2.32| of the ingoing radial (BPT) solution. We have obtained the following 
expansion for the radial coefficient by carrying out small-d asymptotics of where 

are the incidence/transmission coefficients of i?™ (i.e., with the specific normalization used in Sec III which is the 
same normalization as in MST). The expression for the coefficient is given in Eq.(3.12) whereas for we 

used Eq.l 68 ST (which requires Ki, and We obtain: 




sgii7r(t-s+i)p('2£ + l)r( 2 ^ + 2 )r(l - s) 1 
r(£ + i)r(£-s + i)r(£ + s + i) 


1 + 


i Hf. 


l-s 


, iiT \ i (£{£ + 1) + s"^) 

+ H,- 2H_s - ln(2(h) + - - j - ” 2£(£ + 1)' ' 


UJ + 


(6 (£2 + £ - 1) s 2 + £(£ + 1 ) (i 5£2 ^ ^ ^ +Hi- 2H2t - 2H2e+i + ln(2w) - f + 7is) 

2£(£ + l)(2£-l)(2£+l)(2£ + 3) 


(5.4) 


2 (y£(^£ + 1) (^~He-s - He+s - He + 2 iJ_s + ln( 2 a;) - y + 7 b ) - 


ITT 


-(^-iI,_«-i/,+,-i7, + 2i7_, + ln(2w)-y+7Bj" + ^-i/E-^t+s-^f; ■ g 

(8^3 _ 16^2 _ 24£ + 9) s4 g^4 ^g£3 _ 4^2 _ 13^ + 4 - 48^3 + lOR^ + 21£ - 18) 


r(2) 


r( 2 ) 


^ 1 - 


4^2(^+l)(2^-l)2(2f + 3)2 


4(2^- 1)2(2£ + 3)2 


4 £ 2 (£+ 1 ) 2 ( 2 ^- 1 ) 2 ( 2 £ + 3)2 


w 


We note that the leading order of Eq.(5.4) agrees with Eq.3.6.13 |44], which is obtained via an independent method 


based on Page ’sm As a token example of the more simplified form that adopts the expansion for a particular value 
of s, we give Eq.(|5.4|) specifically for s = 0: 


pin 


2-^e5"(^+i)r(2£+l)r(2£ + 2) 1 


r3(£ + i) 


U} 


1+1 


+ 


. 2 -£-igi*,r(«+i)r( 2 £ i)t{ 2£ + 2){6ii{£ + 1 ) - 2 ln( 2 w) + itt + Aje - 1 ) 1 


T3{£ + 1) 


(ln( 2 ia;) — iir — 2 ^e) 


15^2 + IM - 11 


23^e5“(^+i)r(£+i)r(£+|) 

7rr(7+l) y(2f_ i)(2£+1)(2£ + 3) 

( 15^2 (3^(£ 1) _ 2 ^{ 2 £ + 1) - 2 i ;{ 2 £ + 2) + 27^) 2 £ + I 


+ 0ij;{£ + 1 ) - I + 


(2£-l)(2^ +l)(2^ + 3) 4(2£-1)2(2£ + 3)2 

0'il){£ + 1)2 + 3tp{£ + 1 ) + 3^/>(^^(£ + 1 ) - (ln( 2 fw) - in- 27^)2 _ ^ 


s = 0 . 


(5.5) 


Here, the function ijiiz) = r'(z)/r(z) is the digamma function and its nth-derivative. We note that the 

digamma function may be expressed as = He-i ~ 7 b, in terms of the harmonic numbers used above. 
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C. BC strength 


In this section we will provide a small-tr for the BC strength . The solution of the BPT equation has the 
following asymptotics, from Eqs.(3.39) and (3.30) (taking the upper sign), 

R'i 


1 




-e 


00 , 


where 




T^y.ira/ \ _ Ay —i 

(oj) = Aj^uj e 


and A'^ is given in Eq.(3.41|. Comparing with Eq.(2.33l, it follows that 

= _«i?7(r, -uj)A-\ 




T~)U,tra/ \ 

R+ (w) 


(5.6) 

(5.7) 

(5.8) 


Now, from Eq.(3.18) together with the analytic continuation of the irregular hypergeometric [/-functions on the 

i?;;(r,a;). 


complex-w plane (Eq.13.2.41 [46]), we find: 

i?“P(r,a;e2’^0 = e-^™i?“P(r, w) -k 

Using the definition Eq.( 2.36|) of the BC strength, it then follows that 


— e 


q'^(cr) = i 


1 - e 


2-K{uj — iu) 


rtv.tra 

R^ 


i-ia) 


We can therefore calculate the ratio of coefficients on the right hand side via: 


^+,,2.--2.. 
— — 7 —^ ^ 1 




(5.9) 

(5.10) 

(5.11) 


whe re we have made used of Eqs.(5.6) and (3.20). From Eqs.(5.10) and (5.11) and expanding Eq.(3.41) for A'^ and 
Eq.(3.19) for A'^, we find the first hve leading orders for the BC strength: 


q^(cr) = 


ct2«+i 27r(-l)^+«r([-s-f 1) 


M2« 
Hi-s + + 


r([-hs-ki) 


1 + 


16/3 -H 39/2 + ii£_ 17 

+ + 


35^ 


£{e + l){2£-l){2i+l){2i + 3) 2{2£-l){2£+l){2£ + 3) 2/(£-f 1)(2[ - 1)(2/-h 1)(2/-h 3) 


2 ln( 2 cr) — 2je 


cr-l- 


—2Hi-s — 2Hi+s — 


Q{e+£-l)s^+ £{£ + 1) (16/3 -H 39/2 11£ - 17) -k 35 ^ 


£i£+l)i2£-l){2£ + l){2£ + 3) 


+ 47b )( 7b - He-s + ln( 2 cr) - - 


2 “ 2 ) ^ 2 ~ ^ 2 ( 

o(d-)^®+^. 


s (7/2 + 7£ + 3s2 - 6) 


7 / 77 ( 2 ) _ jt{2) , 

2/([-fl)(2/-l)(2/-h3) 


+ 21n2(2d)-- 
6 




(5.12) 


Using Eq.(2.44) to relate BPT’s with RW’s < 7 , it is easy to check that the leading order in Eq.(5.12) agrees with 


Eq.41 [5S| (the ‘ig’ here being ‘if’ in |5S]) for all spins s = 0, —1, —2 (the leading order of RW’s q is independent of 
the spin). 


Again, as a token example of the more simplified form for a particular value of s, we give Eq.(5.12) specifically for 
s = 0 : 


q^ (cr) = q{a) = f 7 r(-l)^a; {2 + 


(16/3-^39/2-^11/-17) 


8/3 -/ 12/2 _ 2 / - 3 
'2 (16/3 -H 39/2 117 _ 17 ) (V;(/-H) - ln(2d)) 8/3 -H 27/2 -h 13/ - 14 


-I- 4ii/'(/ -|- 1) — 4i ln(2cr) ) w— 


(5.13) 


( 2 /- l)( 2 /-h l)(2/-h3) 


( 2 /- l)(2/-M)(2/-h3) 


+ 4(V^(/ + 1) - ln(2d))2 - + 


o(a;3), s = 0 . 
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D. l\W'^ 


We are finally in a position to give an expansion for the main target of this section: the radius-independent quantity 
in the BC integrand Eq.(2.19). From Eqs.(|5.4[) and (5.12) it follows that 


-2f+2 


= M 


2s 


(-l)^+s^22-2<;+2.+ lr3(£ _ g -f S -h 1) r_2f+i 


r2 (£+ I)r 2 ( 2 ^ + 2 )r 2 (i-s) 


8H_s - ^ He.s - 2He+s - 4Hi 


2{l -f 1)(£(8€ -h 7) - 6 )s2 l{£ + + 63) +7)- 23) -f 3s^ 


£{£ + l){2i - l){2i + 1){2£ + 3) 


2{l + 1) {8f + 7l-6)s^ + l{l + 1) (32/3 03^2 _ 33) + 3^4 

/(/ + 1)(2/- l)(2/-hl)(2/-h3) 


- Qcr 


2e+3 


+ o{a^^+% 


(5.14) 


where 

^ _ (-6 (£2 + ^ - 1 ) s 2 + £(£ + 1)(11 _ 15^(£ + 1 )) _ 3^4) ^ 2H_s + ln(2d) + 7 ^) 

— 


4i7_, 2i7_ 


£{£ + l) 


£{£ + l){2£ - l){2£ + l){2£ + 3) 
^,2 


+ 2 + -2He + 4i7_, 


£{£ + l) 


1 


+ 2 ) {He-s + H(.+s + 2,He) — - {He-s + He+s) + 


3H, 


( 2 ) 


H, 


( 2 ) 

l+s 

2 


HP - 4:HP - 


3s^ 


(7^2 + 7£ 6) s 


' 2 '2 ' 2£^{£+l)^{2£-l){2£ + l){2£ + 3) 2£(£-h 1)(2£ - 1)(2£-f 3) 

(128£® -h 768^3 1088^4 -h 120^3 _ 354^2 45^ ^4 _^i2£e _ 2016£3 - 1616£4 -h 1472^3 ^728^2 - 722£ + 59 

4^2(^ + 1)2(2£ _ 1)2(2^ + i) 2(2£ + 3)2 + 4(2^- l)2(2£-f l)2(2^-h3)2 

(256£3 -h 1464£^ -h 2660^® -h 1082£3 - 1361£4 - 926^3 i22£2 + 93^ + is) 52 3^3 ^2 


2^2(^ + 1)2(2£ _ 1)2(2£ + 1)2(2£ + 3)2 


2£(£-hl)(2£-l)(2£-f 3) 3 


(5.15) 


H ere w e give the value of g/|VE|2 specifically for the case s = 0. From Eqs.(5.5) and (5.13), or equivalently from 
Eq.(5.14) with s = 0 (or, equivalently, from [2IQ, we have 


q^{a) q{a) {-lYn {2£ + 1)£\ 


\WTf |1E|2 22^-3 I ((2£-h 1)!!)^ 


5-2«+1 _ ^2^+2 


-32£3 - 63^2 _ 7^ + 23 
2(2£-h3)(2£-h 1)(2£- 1) 


+ 477^ 


(-l)^7r / (2£-f 1)£! 


22^-1 

4f 


((27-hl)!!)^ 


-2^+3 


4(1572 157-11) 


' (ln(2cr) + He- 4:H2e + 1e) 


(27- 1) (27-7 1) (27-7 3) 

(-877.2 ^ g^^ ^ 3^(2) ^ 277^2)) ^ 5127^ + 201673 + 16167^ - 147273 - 112872 + 7227 - 59 


(5.16) 


(27- 1)2(27-7 1)2(27 -7 3)2 


-7 o(tT2^''-3)^ s = 0, 


where 77) is the 7-th harmonic number of order r. 

We note the appearance in (5.14) of a logarithmic behaviour in a at order ct- 2^+3 for small-frequency. It is worth 
pointing out the following ‘curious’ fact. Eveii though the logarithmic behaviour appears al read y at second leading 
order both in the BC strength (see Eq.(5.12)) and in the radial coefficient s7?“’“‘’ (see Eq.(5.4)), there is a delicate 
cancellation between the terms which leads to the logarithmic behaviour for appearing not at second order 


but at third leading order instead. As shown later in Sec |VII[ this implies that the logarithmic behaviour of the Green 
function or a field perturbation will also appear at third - as opposed to second, as one might have expected - leading 
order for late times. 


E. Cases 7 = —s 


As pointed out at the start of this section, the expression in Eq.(5.16) might not be valid for the case s = 7 = 0. 


In this section we present the results obtained by carrying out a specific calculation for this case by setting s = £ = 0 


^ We note a typo of an extra overall ‘-1’ in Eq.8 m- 
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right at the start of the calculation. Although the BPT results for £ = — s = 1, 2 can be obtained directly by putting 
in these values into the general s and £(> 0 ) expressions found above, we also include these cases for completeness 
and because they require an extra step in order to obtain the RW results from the BPT results. 


1. Case s = £ = 0 


We have carried out a specific calculation for the case s = £ = 0 by setting these values right from the start and we 
have obtained: 


= -^+ ln(2cr) - ITT + 7E + ^ + 

IVl (7 Z 

— a (15 tt^ — 6 ln(2CT’) (2 In (2a) — Gin + 6je — 8 ) + 36i'yE'^ — 48i7r + 6(8 — ‘(i')E)lE — 82) + o{a). (5-17) 

36 ^ ^ 


We note that the first two leading orders agree with Eq.(5.5) with £ = 0 but the third leading orders differ slightly in 
the term which does not contain any ‘ln( 2 d)’. 

From Eq.(5.17) we readily obtain (assuming evaluation on the NIA, i.e., ly > 0) 


\W^\^ = = 1 _ (2 (ln( 2 d) + je) + l)d+ 

^ ^72 (ln(2d-) + - 60 (ln(2CT) + 'Je) + + 173^ + o (ct^) . (5.18) 


For the ‘BC strength’ we obtain 


— ‘2 

q^{cr) = q{a) = 27rd + ^ (177r — 127r (ln(2d-) + je)) + 

O 

— {I2ln\2a) + 2 ( 127 £; - 17) ln(2d) + 2jEi6lE - 17) + 20 - tt^) + o (d^) . (5.19) 


We note that if instead of calculating the BC strength direct ly for s = £ = 0 we had used Eq.(5.12) (or Eq.(5.13)) 
with s = £ = 0 we would have obtained the result in Eq.(5.19l but with the number 14 in the place of the number 20 
in the 3rd leading order. 

Finally, from Eqs.(5.18) and (5.19) it follows that 


^ ^ \ = 2707 H---— (l27r^ — 132 (ln( 2 ( 7 ) + ■je) — 85) + o (cr^) , £ = s = 0. 

\W'^\ 3 18 ^ \ \ ) ityj I V / ’ 


(5.20) 


The third leading order differs from Eq.(5.16) in that the ‘-85’ is a ‘-I-59’ in Eq.(5.16); this is a consequence of the 
corresponding discrepancies in both q{(j) and A“. 


2. Cases I = —s = 1,2 


By inserting the values £ = —s = 1 and £ = —s = 2 into Eq.(|5.14|) we respectively obtain 

-a^ -I- o (d®) , £ = —s = 1, (5.21) 


AM 


2 q^{cT) 477 _3 22777 _4 77 (- 1200772 -h 5640 (ln(2d)-£7 B)-£ 10457) _5 


\WT\ 


2 - + 135 


8100 


and 


^ ~ 19177^6 77 (19600772 - 59920 (ln(2d)-£ 7 b) - 348727) 


75 3500 


a"^ + o (ct^) , £ = —s = 2 . 


4410000 


(5.22) 
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VI. RADIAL SOLUTIONS 


In this section we obtain the small-frequency behaviour of the upgoing and ingoing radial solutions valid for arbitrary 
radius r. For this purpose, we use a special trick based on the Barnes integral representation of the hypergeometric 
function. As shown in SecjlTj the radial-dependence of the BC contribution to the GF only comes in through the 
ingoing radial solution, not the upgoing one (see Eqs.(2.21) and (2.37)). Therefore, the radial-dependence at late 


times of the GF itself also only comes in through the ingoing radial solution. For this reason, we apply the mentioned 
Barnes trick to give explicitly the first three leading orders of the ingoing radial solution, whereas we only give the 
leading order (for which the Barnes trick is not necessary) of the upgoing solution. We show, however, how the Barnes 
trick can be used to obtain the behaviour of the ingoing solution up to arbitrary order in the frequency within the 
perturbative regime and we note that it could be similarly applied to the upgoing solution. 

The technique we use for obtaining the small-frequency expansion of the radial solutions can be applied just the 
same to the RW or to the BPT solutions. Besides, from the expansion for the RW solutions one can readily obtain 


the expansion for the BPT solutions via the Chandrasekhar transformation of Sec IIC For this reason, we only give 
the expansions explicitly for one type of solutions: the BPT solutions for the upgoing modes and the RW solutions 
for the ingoing modes. 


A. Upgoing radial solution 


The results in this subsection correspond to the particular normalization of the BPT solutions that is chosen in ST. 
Taking the small-w asymptotics of Eq.4.9 MSTa we obtain 


ROP = 


V2T{2i+l)T{2t + 2) 
2^r(£-h l)F2(£-hs-f l)a;^+«-i 


+ o (w-^-^+2) 


( 6 . 1 ) 


where we have used Eqs.(5.1) and ( |5.3[ ). From Eq.(3.25) for R'^ we then obtain the leading-order asymptotics: 

V2T‘^{2i+l)T{2i + 2 ) (r-l)-"2Ei (-£,-£-s;- 2 £;i) ^ 


Ry = 


2^r2(£ + i)r 2 (£ + s + i)r(£ -s + i) 


/ .^+s+1 


CJ- 




( 6 . 2 ) 


B. Ingoing radial solution 


In this s ubsec tion we consider the ingoing radial solutions of the RW equation. We start with Eqs. (2.16), (3.12), 
dnj ) and ( |3.4a| r 


/f(r,w) = — 


A” 


)VinV°° r(a)r(b) 

2^k——oo 


(6.3) 




^ r(a)r(b) 

Z^fc=-oo 


fe=-oo r(c) 

oo 


— r(c) k— — (x> 


r(a)r( 6 ) f ^ ^ 
^ r(c) 


where 


a = k + v + s + l — ioj, b = —k — v + s — iio, c = 1 — 2iu:. (6-4) 

Let us here carry out a basic comparison of MST’s leading order behaviour with other results in the literature. 
We can use the small-frequency asymptotics for the ‘in’ solution of the Teukolsky equation given by Eqs.3.6.13 and 
3.6.15 [33] (see also |3SlllZj)- Setting a = 0, we obtain: 

sR}t{r, uj) - - l)-"-“V-"+“ 2 Ei(-£ - s, £ - s + 1,1 - s - 2wq 1 - r) = (6.5) 

= ’‘’■“(f - 1 )-®-“V-‘=* 2 Ei(-£ - 2Iw,^ -h 1 - 2iiu, 1 - s - 2wf; 1 - f), (f - l)w < £ -k 1. 


8 

9 


We note that there is a typographical err or o n the right hand side of Eq.(6) Elj: there is a factor missing and the j-s\im in the 

denominator is missing as in Eq.(6.4l but with k ^ j. 

Note that there is a typographical error in Eq.3.6.15 M- a factor r is missing on its right hand side. 
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In order to obtain the RW ingoing solution from BPT’s one, we need to apply the Chandrasekhar transformation as 
per Eq.(2.40). Specifically setting s = 0 in Eq.(6.5), and taking into account of the normalizations as per Eqs.(2.11) 
and (2.32), we obtain 


fe{r, uj) - £+1,1- 2cDi; 1 - f), (f - l)w < £ + 1, s = 0. 


( 6 . 6 ) 


The leading order (w = 0) behaviour of Eq.(6.6) clearly agrees with that from Eq.(6.3). Note that the asymptotics of 
Eq.(6.5) are not amenable to carrying out Fourier frequency-integrations. We now proceed to give a prescription for 
expanding the ingoing RW solution to arbitrary orders in the frequency. 

Now that we know from SeclIVI the behaviour of the coefficients Ofc as a) 


0 our challenge is obtain a suitable 
expansion for the hypergeometric functions in Eq.(6.3). To this end we employ the Barnes integral representation 
(e.g., Eq.15.6.6 


]) which gives 

r(a)m 

r(c) 


2 Fi(a,b,c,z) = ^ J 


r(a + t)r{b + t)r{-t) 
r(c -f t) 


{—zfdt, a, 6 ^ 0, —1, —2,... 


(6.7) 


where the path of integration must be chosen to separate the poles of r(—t) at t = 0,1, 2,... from those of r(a -I- t) 
and r(& -I- t), respectively at t = —a, — a — 1 , —a — 2,... and t = —b, —6 — 1, —6 — 2,... 

For small <D we know from Eqs.(4.1) and (4.2) that 

i/ = £ + 6iy = £ — 1^2^^ -I- O(w^), 

where 1^2 is positive so that the corresponding poles lie at unit intervals left from 

—k — £—s—l + iu} + V2Up' + O(w^), 
k + £— s + iuj — V 2 dP' + 0{d^). 


Note that when a; = 0 we have double poles and it is not possible to find a splitting contour and the representation 
breaks down. 

In Fig. [2 we illustrate the poles and contour Ci of the Barnes integral representation for the term of the MST series 
Eq. ( |6.3[ ) for fc = 1, £ = 2, s = 0, w = —0.1b The proximity of the poles for small uj leads to threading contour Ci 
complicates evaluation of the integral however we may overcome this simply by deforming the contour to the contour 
C 2 lying to the right of all the poles of r(a -I- t) and r(& -|- t) collecting the residues from the poles of r(—£) as we 
do. The contour C 2 can now be taken well away from any poles for example along the line 7/2 + it in Fig. [Band is 
exponentially convergent. 

Letting TV = max(fc -f £—s, —fc — £—s — 1) denote the rightmost pole ofr(a-|-£) and r(&-|-£), the sum over residues 
gives 


N 


Fi{k) = Y. 


n—0 


(— 1 )” r(n -|-fc-|-£-|-s-|-l — iw-l- 5v)T{n — k — £ + s —ioj — Sv) 

r(n -I- 1 — 2iLu) 


{f-iy 


( 6 . 8 ) 


where the coefficients of the polynomial may readily be expanded about w = 0. The corresponding contribution from 
C 2 is given by 


F2{k) = ( 6 . 9 ) 

00 

1 f r{N -I- fc -I- £ -I- s -I- 3/2 — iui + Siy + iy)T{N — k — £ + s + 1/2 — iui — 6iy + iy)T{—N — 1/2 — iy) -.ffj^i/ 2 +iy, 

T{N + 3/2-2ioj + iy) ~ ^ 

— 00 


From the above, we obtain 


feir,uj) = r 


i;s+l 


^(r-1)- 


E OO 

k— — c 


r(a)r(6) 

> “fc r(c) k=- 


^ ak-iFyk)+F2{k)). 


( 6 . 10 ) 


We note that the integral in F 2 is exponentially conv erge nt. All th at one has to do is completely trivial small-a) 
expansions of the various coefficients appearing in Eqs.(6.8l and (6.9). For that, the small-w expansions for v and for 
Qk of the previous section are required. 
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FIG. 1. (Color online). Plot of the contour deformation in the complex t-plane for th e Ba rnes integral representation Eq.( 6.71 
for the hypergeometric function. The straight blue curve is the contour ( Ci ) in Eq.(|6.7[). The dashed red vertical line is the 


contour ( (72 ) in Eq.(6.9l. The circled crosses correspond to the poles of the various T functions in the integrand of Eq.(6.7l: 
the blue ones (at t = 0,1, 2,...) are those of r(— t); the green ones (at t = —a, —a — 1, —a — 2,...) are those of r(a + t); the 
red ones (at t = —b, —b — 1, —6 — 2,...) are those of r(6 + t). This plot is for the case k = 1, i = 2, s = 0, to = —O.li, which 
yields a = 3.9 + 0(a)^), b = —3.1 + 0(lu^) and N = 3. 


C. Cases i = —s = 0,1, 2 


The results of the previous subsection are valid for any values of the spin s (except for, obviously, Eq. 
multipole number |s|. In this subsection we use those results to give the small-tc expansions for the ingoing RW 
solutions specifically for the modes i = —s = 0, 1 and 2, as these are the modes that we will use in the next section 
in order to calculate the Schwarzschild black hole response to a specihc perturbation. 

Specifically, for I! = s = 0 we have, for general r. 


(6.6)) and 


fi{r, -ia) = r + r < - log (r - 1) - 1 + 


^ (r — 1) 2+™sech(7ru) 1 


1 + 2iv 


-|r^ + 8r — 6 + 3 log (r — 1) (log (r — 1) + 2) — 

6 ^ 


_ 3 , _ 3 . 3 . +i_ _ 5,3 + 


(4u^ — 8iv — 3) 


1 + 2iv 


cr^ + o (cr^) . (6.11) 


For £ = —s = 1 we have, for general r. 


fe{r, -ia) = | ^ (3 - 2r + r log (r - 1)) + ^ ^ dv {r - l)2+™r (-iv - |) E (iv +\)\<y+ 

20 — 168f + 144f^ + 16r^ + 3f‘* + 15r log (f — 1) (6 — 4f + r log (r — 1)) + 

1 f°° (f — l)5+™sech(7ru)r (iu + i) , ^ ^ ^ 

dv- -irw-ax—- - (809 — 15r(l + 2iv) — 60(7 + 2iv) log (r — 1) + 


r(*u + |) 


154w 


720 480 240 M d 


2iv + 1 2iv + 3 2iv + 5/1 30 


XX + o o- . 


( 6 . 12 ) 
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For £ = — s = 2 we have, for general r, 
fe{r,-ia) = 

f 1 12 r'^ g ^ 

— < - (6f^ log (f — 1) — 19f^ + 30f^ — 15f + lO) H-z / dv{f — 1)5+™ p (—| — in) F (^ + iv) > d+ 

L 0 J — OO J 

I [840f log (f — 1) (3f^ log (f — 1) — 19f^ + 30f^ — 15f + lO) + 

5040r 

360r® + 1728f® + 40385r^ - 67020f^ + 494707^ - 3977607 + 24860] 

1 f7 — 1)5+™F fin + 4) 

+ / dv- -t: -jgT-^ [1757r(7 — 1)(13 + 2in)sech(7rn)+ 

ooTrr J—QQ T [iv H— 

F (-| - in) r(^ +in) (4207+ 420 log (7- 1) +4201/^ (in + i) - 4201/^ (in + ^) - 1019)] + o (ct^) , (6.13) 


where i/j is the digamma function. 


VII. LATE-TIME TAIL 


In this section we illustrate how one can apply the small-frequency expansions of the previous sections to obtain 
physically-relevant results: the late-time behaviour to high order of the black hole response, at an arbitrary point, to 


a field perturbation of arbitrary integer spin. Before we do that in Sec VIIB we first derive the late-time tail of the 
Green function in the next subsection. 


A. Late-time Tail of the Green Function 


The late-time behaviour of the RW Green function is, via Eq.(2.17), given by t hat of its f-modes sGi- In its turn, 
the late-time behaviour of sG™* is dominated by the BG f-modes sGf^ in Eq.(2.19). Finally, it follows from the 
latter equation that sGf" at late times is given by the small-frequency behaviour of dgCg. We now give more detailed 
expressions for these expansions. 

The radiu s-independent part of th e BG i ntegr and in the BPT case is given in general (except for £ = s = 0) 
in Eq.(5.14), and specifically in Eqs.(5.20)-(5.22) for the lower modes. Let us generically write its small-frequency 
expansion as: 


\WT\- 


= i;ocr^^+^ + i;itT^^+^ + (v2a + V2bina) + o((t^^+^). 


(7.1) 


where the ct- independent constant coefficients vq, vi and V 2 a/b £^re readily readable from the mentioned equations. 
In order to obtain the corresponding RW quantity, we can trivially use Eq.(2.46). We write the small-frequency 
expansion of the proportionality constant in Eq.(2.46l as 


C — Cq + ClfT + C 2 O' + 0(fT ) 


(7.2) 


of the ingoing RW solution as 


where the (T-independent constant coefficients Co/ 1/2 can be read off from from Eqs.(2.47) and (2.481. 

The radius-dependent part of the BG integrand in the RW case is given in general i n F!q .( 6T0| ; all that one h as to 
do is a tr ivial small-a) expansion of the s umm and in Eq.( |6.8| ) and the integrand in Eq.(6.9) (with the use of Eqs.(4.1) 
and (4.31 in general, and specifically Eq.(4.101 for the lower modes). Again, let us write the small-frequency expansion 


fi{r, -ia) = feo[r) + fii{r)a + fi 2 (r)(j^ + 0(cr'^), 


(7.3) 


where the (T-independent (but radius-dependent) coefficients 0/1/2 can be readily obtained in the manner just 
indicated. 

T he an alytic small-frequency expansions of the radius-independent (Eq.(7.1| times Eq.(7.2|) and radius-dependent 
(Eq.(7.3) evaluated at r times the same expression but evaluated at r') parts of the BG integrand in the RW case are 
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then to be put together in Eq.( |2.21[ ) and integrated as per Eq.( |2.19[ ). The result, for late times is, straight-forwardly, 

(7.4) 


■sGf^ — i CoVofooifTf') + (2^* + 3) 




+ {CqVi +CiVq) fQo{r,r') t ^ + 


2(^ + 2)(2^ + 3) 


Co^^o/ol 2 (?', r') + (vqCi +CoVi) foi{r,r')+ 


{CoV 2 b ('0(2^ + 5) -Int) +CQV 2 a +U 0 C 2 +C 1 V 1 ) /oo(r, r') 




where 


foo{r,r') = fm{r)fmir'), foi{r,r') = fm{r)fii{r') + feo{r')fii{r), 
foi 2 {r,r') = feo{r)fe 2 {r') + fm{r')fe 2 {r) + fn{r)fa{r'). 


(7.5) 


Here we are using the dimensionless time t = t/(2M). The leading-order behaviou r at late times is, therefore, 
sGf^ ~ for all integer spins. We also note that the logarithmic behaviour in Eq.(5.14) for /\W'^\^ for small 


G?^ as t ® I nf. T he late-time behaviour of th e 
in Eq.(2.17) by the expansion in Eq.(7.4). 


frequencies led to the appearance of a logarithmic behaviour in 
GF of the 4-dimensional RW equation is then obtained by replacing 

A similar analysis can be done for the BPT GF, but using the corresponding equations instead: (2.37) and (2.35). 
The radius-independent part of the integrand, we have already given for the BPT case. The radius- 

dependent part can be obtained from the expansions Eq.(7.3) of the ingoing RW solutions via the Ghandrasekhar 


transformation of Eq.(2.39). The late-time behaviour of the f-modes of the BPT GF is of the same leading order as 
that in the RW case (i.e., f“^^“^) and the logarithmic term also generally appears at the same order as in the RW 
case (i.e., Inf). 


B. Late-time Tail of an Initial Perturbation 


We shall now give a particular application of our late-time results for the GF. Let us here consider an initial field 
perturbation given by = 0) for the Gmultipole of the field and = 0) for the 

Gmultipole of the time derivative of the field. Then the response of a Schwarzschild black hole is given by 


ue{r»,t) = 


dr' 


,Gr (r, r'; tWir'J + «-(r:) ^^Gf (r, r'; f) 


(7.6) 


where GY^*{r,r';t) is the Gmultipole of the retarded GF of the wave equation satisfied by the field. In this section 
we will take the RW equation as th e wave equation. The late-time asymptotics of the perturbation ug are given by 
replacing ^G)®* by sGf^ in Eq.(7.6) and approximating sGf^ by performing a small-frequency expansion of SsG^ in 

Eq.lirM. 


As in m let us consider the following initial perturbation for general spirp^ 


^cr \ 1 ( ~ (^* “ ^ 0 ) 


2M2 




(7.7) 


with xq = r*(10M). 

In Figs|2]^we plot the time evolution via the RW equation of a spin-s field for various multipole numbers i using 
the initial data Eq.(7.7), similarly to m where we only presented the case s = 0, i = 1. We calculate the time 
evolution numerically up to T = 3000M, where T = t — 2xo, and we compare it with late-time asymptotics. The 
numerical solution is obtained using the code in |49) for the (1 -I- l)-dimensional differential equation which results 
from the 4-d RW equation after factorizing out the angle dependence of the solution via scalar spherical harmonics. 
We obtain the late-time asymptotics, up to four leadi ng or ders, from the results of the previous section. In order to 
obtain the small-frequency expansion of SgG^ via Eq.(2.21), we need the small-frequency expansion of g/|IFp and of 


As an alternative to the initial data of Eq. ( ( 7 ^ we could consider the following initial data (as used in m) : zero for the field and a 
Gaussian distribution of small width for the time-derivative of the field. In this case, the perturbation response is an approximation 
to the Green function. However, we found that such approximation is not as good as 113 for assessing the validity of the late-time 
asymptotics to the higher-than-leading orders as intended in this paper. 
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the ‘in’ radial solution fi. We now give the specific expansions of the perturbation response for different spins s and 
i modes. 


We found the late-time asymptotics for the case s = 0 and £ = 0 using Eq.(6.11) for th e radial solution and Eq.(5.20) 
for g/lVEp ; for £ = 1, 2 we merely inserted the corresponding value of £ into Eq.(5.16). The results are the following 


(the approximation sign is due to the fact that we have rounded up the coefficients to seven significant figures): 


U£ 


BC 


20.02467 269.1100 -881.0855 In (T)-t 7561.913 


r3 




2^5 


+ o{f-^), s = £ = 0, 


(7.8) 


gc 1144.318 52385.24 43484.10 In (T) - 1746947 


J'b 






-7 


), s = 0, £ = 1, 


(7.9) 


HP 42634.22 3114525 -1796322 In (T)-t 146090980 , _ 

-+-y- +o{T-^), s = 0,£ = 2. 


2^7 


2^8 


2^9 


(7.10) 


Eigsj2]^show excellent agreement at late times between the numerical solution of the RW equation for s = 0 and 
= 0,1, 2 and the late-time asymptotics of Eqs.( 7.8)-(7.101. 


Mu, 



FIG. 2. Scalar field Mue at r* = xg as a function of time T/M for the initial data Eq.(7.7l with xg = r*(10M) and s = £ = 0 


(i.e., similar to Fig.l in m and with the same initial data bnt here it is for s = £ = 0). Red curve: numerical solution. Dashed 
black curve (overlapping wit h red curve): late-time expression given by Eq.(7.8l (which comes from Eqs.(5.20l together with 
the small-CT results in SecfVI for the radial function). Green curve: numerical solution minus the leading order (i.e., 0{T~^)) in 


the late-time expression. Solid blue curve: numerical solution minus the two leading orders (i.e., 0{T~^) and 0{T~'^)) in the 
late-time expression. Solid cyan curve: numerical solution minus the four leading orders (i.e., 0{T~^), 0{T~^), 0(T“® ln(T)) 
and 0{{T)~^)) in the late-time expression. 


The late-time asymptotics for the RW cases |s| = £ = 1 and 2 are obtained similarly to those above for s = 0: we 
used Eqs.(6.12) and (6.13) for the small-frequency expansions of the radial soluti on an d Eqs.( 5.21 ) and (5.22) for that 
of /\W-^ which is then simply converted to RW version with the use of Eqs.(2.44) and (B4|. The final results for 
the late-time perturbation are the following: 




ia70216-028M.051n(f) 


and 


Ui 


BC 


60748.13 4318418 (198106997 - 3466693 In (T) 


2^7 


2^8 


2^9 


+ o(r-®), |s|=£ = 2. 


(7.12) 
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MU/ 



FIG. 3. Similar to Figj^but here it is for s = 0, ^ = 1 (the orders subtracted from the numerical solution are the cor resp onding 
leading order, next-to-leading order and four leading orders for this case). The late-time expression is given by Eq.(7.9|. 


Mu/ 



FIG. 4. Similar to Figj^but here it is for s = 0, ^ = 2 (the orders subtracted from the numerical solution are the corresponding 
leading order, next-to-leading order and four leading orders for this case). 


Figsshow excellent agreement at late times between the numerical solution of the RW equation for |s| = £ = 1, 2 
and the late-time asymptotics of Eqs.(7.11)-(7.12). 


The above examples illustrate the appearance of a higher-order logarithmic behaviour of a perturbation response, 
appearing after Price’s well-known power-law tail Hill- We note that, as already pointed out in Sec|VD[ this 


logarithmic behaviour appears at third leading order at late times. It appears at third leading order, not at second 
leading order as one might have naively expected from the fact that the logarithmic behaviour at small frequency 
appears at the second leading order both in the BC strength and in the Wronskian. This is due to there being 
delicate cancellations between certain terms of and the Wronskian, which we have been able to derive using the 
precise values of the coefficients in the expansions. 
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MU/ 



FIG. 5. Similar to Fig|^but here it is the RW solution for |s| = £ = 1 (the orders subtracted from the numerical solution are 
the corresponding leading order, next-to-leading order and fonr leading orders for this case). 


Mu; 



FIG. 6. Similar to Fig|^but here it is the RW solution for |s| = £ = 2 (the orders subtracted from the numerical solution are 
the corresponding leading order, next-to-leading order and four leading orders for this case). 


VIII. CONCLUSIONS 

In this paper we have developed the MST formalism for the solutions of the radial Regge-Wheeler and Bardeen- 
Press-Teukolsky equations, which are obeyed by linear field perturbations of a Schwarzschild black hole space-time. 
We have derived, for the first time, the MST formalism for the solutions of the RW equation for spin-1 as well as the 
MST expressions for the branch-cut-relevant quantities for general spin. We have given explicit expansions for small 
frequency up to the first four leading orders for the various MST quantities for general spin. In principle, the MST 
series could be expanded to arbitrarily large order in the frequency. The main difficulty in achieving that is the fact 
that the small-frequency expansion of the renormalized angular momentum parameter iz that we currently use does 
not reproduce the numerically-observed behaviour of a sudden appearance of an imaginary part in v as the frequency 
is increased from 0 to larger real values. 

We have used our small-frequency expansions in order to obtain the late-time behaviour for arbitrary radius of 
spin-field perturbations of a Schwarzschild black hole up to the first four leading orders (for a specihc multipole-£). 











36 


Our results explicitly reveal a new logarithmic behaviour at third order for late times as Int. We note that the 

appearance of a logarithmic behaviour was already predicted by other works (see, e.g., [26115011^ 1. However, to the 
best of our knowledge, the order at which the logarithmic behaviour appears was not correctly predicted anywhere 
(nor was the calculation of the coefficients for general radius carried out). As we noted at the end of Sec V D there 
is a delicate cancellation between different logarithmic terms which averts the appearance of a logarithmic behaviour 
at a lower order; this delicate cancellation is probably hard to predict unless an exact and detailed analysis is carried 
out such as the one in this paper (see also |53j . where they find numerical evidence that the logarithmic behaviour 
does not appear at first nor second leading orders). 

We succinctly presented in the Letter m the final results that we have derived in this paper. We already used 
some of these results in the calculation of the scalar self-force in Schwarzschild space-time carried out in [TS]. The 
natural extension of our results is the explicit calculation of the small-frequency expansions of the MST quantities 
for the Teukolsky equation in Kerr space-time, with the corresponding late-time analysis of perturbations of a Kerr 
black hole. We expect to present the Kerr analysis in the near future, together with its application to a self-force 
calculation. 

Finally, while the calculation of the quasinormal modes has been applied to the modelling of radiation after the 
inspiral of two black holes via a matching to a numerical relativity solution, to the best of our knowledge, the branch 
cut has never been taken into account for such purposes. We would expect that the inclusion in the modelling of the 
branch cut, together with the quasinormal modes, would help match the analytical solution from perturbation theory 
with the numerical relativity one. 
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Appendix A: Validation of BC results 


In this section we compare the small-d results presented in this paper with an independent method which we 
presented in [27]. The latter method is naturally-adapted to a ‘mid-frequency’ regime (which, although valid for all 
frequencies, is not practical to use in the asymptotic small- or large-frequency regimes). In this appendix we illustrate 
with plots that there is a region of overlap between the results obtained with the small-d method presented here 
and those obtained with the mid-frequency method of m- In |?7p] we already presented plots of quantities focused 
mainly on the s = £ = 2 case. Therefore, for variety, here we will focus on the s = 0 and £ = 0,1 cases. 


As explained in [27], the ingoing radial function fe{r,uj) possesses simple poles along the NIA (which are ultimately 
irrelevant, as they are cancelled out in the Fourier modes of the GF by the corresponding poles in the Wronskian). It 
is therefore useful to define the following radial function: fi(r,u}) = sin(—27riw)/^(r, w). 

In Figsl^ 1^ l^and 10 we respectively plot the radial function /^, the BC ‘strength’ g(cr), the Wronskian W and the 
BC mode SsG^r, r', —icr), all as functions of the frequency a along the NIA. 


Appendix B: Radial Coefficients 


In this appendix we use the Chandrasekhar transformation Eq.(2.38) in order to relate the BPT and RW radial 
coefficients and Wronskians. Let us first add higher orders to the asymptotics in Eq.(2.111 of the RW solution: 


(1-I-a+A-I-/3+A^) e r*-)■-oo, 

Af (l -k ^ -k -k , r, ^ -koo. 


(Bl) 


We note that there are two typos in Eq. 3.6JW|: there should be an overall factor ‘—1’ in the expression for A/i^(r, i/), using the notation 
of that paper, and, in the expression for it should be an+ instead of dn-- The minus sign typo carried over to Fig.10 m, which 

is therefore a plot of Ag£\ instead of ^—ie'^‘^*Ag£ as stated there. 
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--/tLir.r 

e fj 



FIG. 7. The radial function /^(r, ai) at r = lOM as a function of a for s = 0, £ = 1. We calculated the straight 

blue curve using the so-called Jaffe series (see, e.g., [Ill 127] '). and the dashed red curve using the small-frequency asymptotics 
presented in this paper up to 0(a^^) (specifically, using Eq.l37 [TH as in this case we did not need to obtain the coefficients of 
the expansion analytically via the use of the Barnes integral method). 


q q 



FIG. 8. BG ‘strength’ q{a) of Eq.(2.20l (or, equivalently for s = 0, q^(cr) of Eq.(2.36l) as a function of a. In the straight 
blue curves we calculated q{a) using the mid-frequency method of [57] (with the radial functions calculated at r = 5M). In 
the dashed red curves we calculated q{a) using the small-frequency asymptotics presented in this paper up to (a) For 

s = 0, £ = 0. (b) For s = 0, £ = 1. 


where a+, /3+, cToo and /3oo are coefficients to be determined. By imposing that these asymptotics satisfy the RW 
equation, we find 


Q^OO 


a+ 


iX ^ A(2 — A) -b — I)a; 

= ~ 8w2 

I (A -b I - s^) 4 _ s4 _ ;^2 y(2 - 16zMa;) -b (2A -b 201 Mw - 3) - 20zMa; 

4M2(4Ma7Tf) ’ ~ 


(B2) 


We note that Ofooj/3oo,/3+ are all real-valued when uj is purely imaginary. 


Using the transformation Eq.(2.38l it is straight-forward to find the following coefficients of BPT solutions: 
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W[f,,g,]*M W[?,,g,]*M 




FIG. 9. The Wronskian W ^fe , as a function of a for s = 0, ^ = 1. The straight-blue and dotted-red curves respectively 

correspond to the real and imaginary parts of the Wronskian evaluated using the mid-frequency method of m with the radial 
functions calculated at at r = 2.8M; the convergence of the series used in the mid-frequency method becomes slower as the 
frequency becomes smaller. The dashed blue and red curves respectively correspond to the real and imaginary parts of the 
Wronskian evaluated using the small-frequency asymptotics presented in this paper up to In figure (b) we plot the 

mid-frequency to larger values of the frequency so as to show its general form; the green-dashed curve here corresponds to the 
absolute value of the Wronskian. Cf. Fig. 13 [27] for s = £ = 2. 


dsGi/M 



FIG. 10. BG mode SsGe (r, r', —ia) of Eq.(2.211 as a function of a at r = r(r, = 0.4M), r' = r'(r, = 0.2M) for s = 0 and 
1=1. In the straight blue curve we calculated tbe BC mode using the mid-frequency method of m and in the dashed red 
curve we calculated it using the small-frequency asymptotics presented in this paper up to 


It then immediately follows, using Eqs.(2.13l and (2.34), that 


(2M, 

I Q^cxd 

= w { 

M/3oo 

I a+ -b 2M^p+ ’ 



Note that FF^ is, by definition, independent of the chosen normalization for and 
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